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Generic higher character Lifshitz critical behaviors are described using 
field theory and ex,- expansion renormalization group methods. These criti- 
cal behaviors describe systems with arbitrary competing interactions. We de- 
rive the scaling relations and the critical exponents at the two-loop level for 
anisotropic and isotropic points of arbitrary higher character. The framework 
is illustrated for the N -vector c/> 4 model describing a d-dimensional system. 
The anisotropic behaviors are derived in terms of many independent renor- 
malization group transformations, each one characterized by independent cor- 
relation lengths. The isotropic behaviors can be understood using only one 
renormalization group transformation. Feynman diagrams are solved for the 
anisotropic behaviors using a new dimensional regularization associated to a 
generalized orthogonal approximation. The isotropic diagrams are treated us- 
ing this approximation as well as with a new exact technique to compute the 
integrals. The entire procedure leads to the analytical solution of generic loop 
order integrals with arbitrary external momenta. The property of universal- 
ity class reduction is also satisfied when the competing interactions are turned 
off. We show how the results presented here reduce to the usual m-fold Lifshitz 
critical behaviors for both isotropic and anisotropic criticalities. 
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I. INTRODUCTION 



Field theoretic renormalization group techniques are invaluable tools for studying usual 
critical phenomena as well as the critical behavior associated to the physics of systems 
presenting arbitrary short range competing interactions. The universality classes of the 
ordinary critical behavior are characterized by the space dimension of the system d and the 
number of components of the (field) order parameter N [1,2]. Competing systems, on the 
other hand, possess different types of space directions known as competition axes. 

The simplest type of competition directions can be most easily visualized using the 
terminology of magnetic systems via a generalized Ising model. One permits exchange fer- 
romagnetic couplings between nearest neighbors (J\ > 0) and antiferromagnetic interactions 
between second neighbors (J 2 < 0) occurring along m-2 dimensions. Whenever mi < d the 
system presents a (usual) second character anisotropic Lifshitz critical behavior whose uni- 
versality classes are characterized by (N, d, m 2 ), whereas the isotropic behavior characterized 
by d = m 2 close to 8 was formerly described at the same time [3]. The phenomenological 
model corresponding to a uniaxial anisotropy (m 2 = 1) in a cubic lattice is known as 
ANNNI model [4]. In the critical region, this sort of system is characterized by a disor- 
dered, a uniformly ordered and a modulated ordered phase which meet in a uniaxial Lifshitz 
multicritical point, where the ratio ^ is fixed at the corresponding Lifshitz temperature 
Tl. High-precision numerical Monte Carlo simulations were carried out for the critical ex- 
ponents of this model [5] and checked using two different two-loop analytical calculations 
[6-8]. From the renormalization group perspective there is an important difference between 
these Lifshitz critical behaviors. The anisotropic behaviors have two independent correlation 
lengths, perpendicular to the competing axes as well as ^4 parallel to the m 2 competing 
axes. The isotropic behavior has only one correlation length £1,4. 1 

If we go on to include ferromagnetic couplings up to the third neighbors ( J3 > 0) along 
a single axis, the system will present a uniaxial third character Lifshitz point whenever 
^ and ^ take certain fixed values at the corresponding Lifshitz temperature [11]. When 
this sort of competition takes place along m 3 spatial directions, the system presents a m 3 - 
fold third character Lifshitz point. On the other hand, if simultaneous and independent 
competing interactions take place between second neighbors along m 2 space directions and 
third neighbors along 777,3 space dimensions, the system presents a generic third character 
7773-fold Lifshitz critical behavior. The generic third character universality classes are defined 
by (N, d, 777, 2 , m 3 ), thus describing a wider sort of critical behavior when compared with the 
thrid character universality classes (N,d,m 3 ). 

This idea can be extended in order to define the m^-fold Lifshitz point of character L, 
when further alternate couplings are permitted up to the Lth neighbors along mi directions, 
provided the ratios ^j^- jj± take especial values at the associated Lifshitz temperature 
[12-14]. However, the most general anisotropic situation is to consider several types of 
competing axes occurring simultaneously in the system such that second neighbors interact 
along m 2 space directions, 777 3 directions couple third neighbors, etc., up to the interactions 
of L neighbors along m L dimensions, with all competing axes perpendicular among each 



For an alternative field-theoretic approach for m-axial Lifshitz points, see [9,10]. 
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other. In that case, the corresponding critical behavior is called a generic Lth character 
Lifshitz critical behavior [15]. 

In this work we shall undertake an exploration of the field theoretical renormalization 
group structure of the most general competing system using e^-expansion techniques for 
anisotropic and isotropic higher character Lifshitz critical behaviors. A rather brief descrip- 
tion of this structure was set forth in a previous letter [15] where it was first described; here 
we shall present the details and extend the formalism in order to incorporate the exact two- 
loop calculation for arbitrary isotropic higher character criticalities. Renormalization group 
(RG) arguments are constructed in order to find out the scaling relations for the anisotropic 
as well as the isotropic critical behaviors. The discussion parallels that for the usual second 
character Lifshitz points [8]. 

The arbitrary competing exchange coupling Ising model (CECI model) is the lattice 
model associated to this new critical behavior. It has a general renormalization group 
structure which contains many independent length scales, and its construction can be utilized 
for both anisotropic and isotropic cases. In the anisotropic cases, the system has only nearest 
neighbor interactions along {d — m 2 — ... — mi), second neighbors competing interactions 
along m 2 dimensions, and so on, up to Lth neighbors competing interactions along tjil spatial 
directions. The distinct competing axes originate several types of independent correlation 
lengths, namely £i for directions parallel to the (d— m 2 — m L )-dimensional noncompeting 
subspace, £ 2 for directions parallel to the m 2 -dimensional competing subspace, etc., and 
characterizing the m^-dimensional subspace. The simplest representative of the CECI model 
is better understood with the help of Fig.l, which is the particular case m 2 = rris = 1, 
?7i4 = ... = rriL = 0. There are two competing subspaces and three types of correlation 
lengths which define three independent renormalization group transformations. It defines a 
particular generic third character anisotropic Lifshitz critical behavior. 

In the phase diagram of the ANNNI model, the parameters which are varied are the 
temperature T and p — ^ which take a particular value at the uniaxial second character 
Lifshitz multicritical point as depicted in Fig. 2. It is a particular case of the CECI model 
whenever m 2 = 1, with m 3 = ... = rriL = 0. Although the ANNNI model has applications 
in several real physical systems (see for example [4]), the prototype of second character 
Lifshitz points in magnetic materials is manganese phosphide (MnP). Experimental as well 
as theoretical investigations have determined that MnP presents a pure uniaxial Lifshitz 
point (m 2 = 1, d = 3, N = 1) [16,17]. 

When adding further competing interactions to the ANNNI model, the number of pa- 
rameters in the phase diagram increases [12]. For instance, in the phase diagram of the 
model including uniaxial competing interactions up to third neighbors the parameters to be 
varied are the temperature T, p\ — ^ and p 2 = j^. One can locate the third character 
Lifshitz point by looking at the projection of the phase diagram in the plane (pi,p 2 ), as 
was demonstrated using numerical means [12]. In the example of the CECI model displayed 
in Fig.l, let the competing exchange be completely independent along the different com- 
peting axes. In that case, the phase diagram can be described by T, p z = ^ p lv = ^ 

P2y = A useful two-dimensional representation can be obtained by separating the phase 
diagrams in two parts. The diagram (T,p z ) characterizing the second character behavior 
(with pi y ,p 2y fixed) and the diagram (pi y ,P2 y ) (with T, p z fixed) corresponding to the third 
character behavior can be ploted independently. The superposition of the two diagrams at 
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the generic third character Lifshitz point is indicated in Fig.3. As a consequence, there is 
a uniformly ordered phase and two modulated phases called Helical and Helical^ in Fig.3 
which meet at the uniaxial generic third character anisotropic Lifshitz point. Now there are 
two first order lines separating the ferromagnetic- Helical 2 and Helical 2 — Helical^ phases. 
Analogously, when there are arbitrary independent types of competing axes, we can consider 
several independent phase diagrams and each two-dimensional projection of them. The su- 
perposition of them in one two-dimensional diagram gives origin to a situation that resembles 
that illustrated in Fig.3 . Instead, there are several modulated phases and one uniformly 
ordered phase which meet at the generich Lth character anisotropic Lifshitz point. Each 
competing subspace has its own characteristic modulated phase along with its own indepen- 
dent correlation length. Although these objects go critical simultaneously at the Lifshitz 
critical temperature, they define independent renormalization group transformations in each 
different subspace. 

Therefore, we find multiscale scaling laws as a consequence of this renormalization group 
flow independence in parameter space. This implies that we find several independent cou- 
pling constants, each one depending on a definite momenta scale characterizing the particular 
competition axes under consideration. Nevertheless, all coupling constants flow to the same 
fixed point. The universality classes of this system are characterized by the parameters 
(N, d, m 2 , rn L ), therefore generalizing the usual Lifshitz behavior. It is important to men- 
tion that when we turn off all the competing interactions between third and more distant 
neighbors, the universality classes of the generic higher character Lifshitz point turn out to 
reduce to that associated to the second character behavior (TV, d, 777.2). Notice that these 
anisotropic behaviors generalize previous lattice models with competing interactions [18] as 
it includes all types of competing axes. The isotropic critical behaviors d = m n have a 
distinct feature in which there is only one type of correlation length £ 4n . Their universality 
classes are characterized by (N, d, n) where n is the number of neighbors coupled through 
competing interactions. 

In addition, we compute the critical exponents at least at 0(e 2 L ) using dimensional regu- 
larization to resolving the Feynman diagrams and normalization conditions (and) or minimal 
subtraction as the renormalization procedures. The computation is realized in momentum 
space. For the anisotropic cases, the Feynman diagrams are performed with an approxima- 
tion which is the most general one consistent with the homogeneity of these integrals in the 
external momenta scales. The isotropic situations are treated using this approximation as 
well, but we also present the exact calculation at the same loop order and make a comparison 
with the above mentioned approximation. 

We present the functional integral representation of the model in terms of a A0 4 setting 
and define the normalization conditions for this higher character Lifshitz critical behavior 
in section II. We show that many sets of normalization conditions, each one corresponding 
to a specific type of competition axes, are convenient to have a satisfactory description of 
the problem in its maximal generality. 

In Section III we present the renormalization group analysis for the anisotropic critical 
behaviors. We construct the several renormalization functions appropriate to each competing 
subspace and study their flow with the various renormalization group transformations. We 
find the proper scaling relations to each competition subspace. 

Section IV discusses the renormalization group treatment for the various isotropic be- 
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haviors. We obtain the scaling relations and show that they reduce to the usual A</> 4 case 
when the interactions beyond the first neighbors are switched off. 

Section V is an exposition of the calculation of the critical exponents for the anisotropic 
cases using the scaling relations obtained in section III. We describe two different ways of cal- 
culating the critical exponents either in normalization conditions or in minimal subtraction. 
We discuss the limit L — > oo and some of its implications. We point out the analogy of the 
Lifshitz critical region with an effective field theory which arises from a recent cosmological 
model including modifications of gravity in the long distance limit [19]. 

Sections VI and VII are an in-depth analysis of the critical exponents for the isotropic 
cases. The critical exponents for the isotropic cases are calculated using the orthogonal 
approximation and the scaling relations in section VI. The exact calculation of the critical 
exponents and the comparison with those obtained from the orthogonal approximation are 
carried out in section VII. 

The particular case corresponding to the second character isotropic Lifshitz critical be- 
havior is discussed explicitly in Section VIII. The resulting critical exponents are shown to 
generalize those obtained previously [3]. We discuss their relationship with those coming 
from the orthogonal approximation. 

Section IX concludes this paper, a discussion of the ideas is summarized and some possible 
applications will be proposed. We calculate the Feynman integrals in the appendices. We 
describe in detail the generalized orthogonal approximation for the calculation of higher 
loop integrals of the anisotropic behaviors in Appendix A. It will be shown there that one 
can obtain the answer in a simple analytical form for arbitrary external momenta scales. 
The property of homogeneity of these integrals along arbitrary external momenta scales 
is preserved. Then, we use the same approximation to compute diagrams for the isotropic 
behaviors in Appendix B. In addition, we perform the exact calculation for arbitrary isotropic 
cases in Appendix C. We also analyse the simple particular case associated to the usual 
second character isotropic behavior. 

II. FIELD THEORY AND NORMALIZATION CONDITIONS FOR HIGHER 

CHARACTER LIFSHITZ POINTS 

The field theoretical representation can be expressed in terms of a modified A0 4 field 
theory presenting arbitrary higher derivative terms due to the effect of competition along 
the different kinds of competing axes. The type of competing axes are defined by the number 
of neighbors that interact among each other via exchange competing couplings. Let m n be 
the number of space directions whose competing interactions extend to the nth neighbor. 
Thus, the m n -dimensional competition subspace will be represented in the Lagrangian with 
even powers (up to the 2nth) of the gradient acting on the order parameter scalar field. 
Thus, the effect of the competition resides in the higher derivatives of the field. 

The corresponding bare Lagrangian density can be written as [15] : 

^=^v (d _ EL2m „)0oi 2 +x:Yiv™^oi 2 a) 

n=2 
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L j L— 1 rt— 1 j 

n=2 n=3 n'=2 

+ ^000 + ^ A o0o- 

At the Lifshitz point, the fixed ratios among the exchange couplings explained above trans- 
late into this field-theoretic version though the conditions So n = r nn i = 0. All even momen- 
tum powers up to 2L become relevant in the free propagator [20]. This condition simplifies 
the treatment of the system since it allows the decoupling of the several competing subspaces 
of Feynman integrals in momentum space. It indeed makes possible to solve these diagrams 
to any desired order in a perturbative approach. Within the loop order chosen, we can set a 
perturbative regime with maximal generality as far as critical behavior of competing systems 
are concerned. So we need the small loop parameter, which is intimately connected to the 
critical dimension of the theory. 

It is instructive to find the critical dimension of this field theory through the use of 
the Ginzburg criterion [21-23]. From the perspective of magnetic systems in the above 
Lagrangian density, to — toL + (T — T^) measures the temperature difference from the critical 
temperature T^. In the mean-field approximation the inverse susceptibility is proportional to 
(T — Tl), but has no longer this behavior when fluctuations get bigger due to the closenesss 
of the critical point and the mean-field argument breaks down. This immediately leads to 
the critical dimension d c = 4 + Yln=2 ^~~^ m n- Above this critical dimension the mean- 
field behavior dominates the system. Consequently, the small parameter for a consistent 
perturbative expansion is = 4 + X^n=2 ^T^ m « ~~ ^- The same type of argument can be 
constructed for the isotropic behaviors. When d = m n , the critical dimension is d c = An— e^. 

The renormalized theory can be defined starting from the bare Lagrangian (1). The ap- 
proach we shall follow here is completely analogous to that exposed in [8] (see also [1]) and 
the reader is invited to consult that reference in order to be familiarized with the notation 
employed. The renormalization functions are determined in terms of the renormalized re- 
duced temperature and order parameter (magnetization in the context of magnetic systems) 

-1 

as t = ZT^to, M = Z£ and depend on Feynman graphs. When the theory is renormal- 
ized at the critical temperature (t — 0), a nonvanishing external momenta must be used 
to define the renormalized theory. Consequently, the renormalization constants at the criti- 
cal temperatute T L depend on the external momenta scales involved in the renormalization 
algorithm. 

Let us analyse the anisotropic behaviors. The Feynman integrals depend on various 
external momenta scales, namely that characterizing the (d — 777,2 — ••• — tul) -dimensional 
noncompeting subspace, a momentum scale associated to the 7772 space directions, etc., up 
to the momentum scale corresponding to the mi competing axes. Thus, it is appropriate 
to define L sets of normalization conditions in order to compute the critical exponents 
associated to correlations either perpendicular to or along the several types of competition 
axes. We define K\ to be the external momenta scale asssociated to the (d — m 2 — ... — Tri- 
dimensional noncompeting directions. If we define the noncompeting directions to be along 
the mx-dimensional subspace, where m\ = d — m 2 — ... — m^, we can unify the language by 
stating that n n is the typical external momenta scale characterizing the m n competing axes 

(77=1,...,L). 
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We now turn our attention to the definition of the symmetry points (SP). In case we 
wish to evaluate the critical exponents along the jth type of competing axes, we set n n = 
for n 7^ j keeping, however, Kj ^ 0. The proper normalization conditions to evaluating 
exponents along the various competition axes can be defined as follows. If k'^ is the external 
momenta along the competition axes associated to a generic 1PI vertex part, the external 

2 

momenta along the nth type of competing directions are chosen as k'^.k'-^ = ^f(ASij — 1). 
This implies that (k'^ + k'-,-,) 2 = K 2 n for % ^ j. The momentum scale of the two-point 
function is defined by k',\ = k 2 1 = 1. The set of renormalized 1PI vertex parts is given by: 



r 2(„)(0^n) = 0, (2a) 
ok/? («) n 

T R(n)( k i(n),9n)\sP n = 9n, (2c) 

^R(n)(K(n)iK(n)ik',g n )\sP n = ^ ( 2 d) 

rS?)(%)^n)L;- = C = o. (2e) 



These L systems of normalization conditions seem to provide L renormalized coupling 
constants. The origin of this overcounting is a consequence of the L independent flow in the 
renormalization momenta scales k\,...,kl- The analysis works with L coupling constants, 
namely g n = u n (K 2 ?)~ (and X n = uq^k^ 1 )^) characterizing the flow along the momenta 
components parallel to each m n -dimensional competing subspace. This is really a disguise 
since the situation becomes simpler at the fixed point: the many couplings will flow to the 
same fixed point, at two-loop level, giving a clear indication that this property is kept in 
higher-loop calculations. The explicit demonstration of this fact will be tackled in Sec. V. 2 

The normalization conditions for the isotropic case (m n = d near An) can be defined 
in a close analogy to its second character isotropic particular case [8]. If k\ is the external 
momenta along the m n competition axes, the external momenta along the An directions are 

2n 

chosen as k\.k^ = ^f-(4<5y — 1). This implies that {k[ + k'^) 2 = for % ^ j. The momentum 
scale of the two-point function is fixed through k' 2n = k 2 ™ = 1. Then we have the same 
normalization conditions Eqs.(2), but now there is solely one type of external momenta scale. 
The others are absent in this situation as an effect of the Lifshitz condition 8o n = T nn t = 0. 

We can express all the renormalization functions and bare coupling constants in terms 
of the dimensionless couplings in a unified perspective for both anisotropic and isotropic 



2 Let us mention briefly another type of anisotropic critical behavior whose theoretical existence 
is granted from the CECI model structure. In magnetic systems, the absence of the ferromagnetic 
phase leads to a structure where several modulated phases are mixed among each other. This 
situation is transliterated in the condition d = mi + ... +raj, In that case, the normalization 
conditions can be defined without making reference to the noncompeting subspace. In general 
the critical dimension will increase and the lowest character modulated phase will play the role of 
the former ferromagnetic ordered phase. We shall not delve any further into this situation in this 
paper, but leave the analysis for future work. 
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behaviors. The subscript n = 1,2,3, ...,L labels the different external momenta scales be- 
longing to the general Lifshitz critical behavior, as defined above for the anisotropic and 
isotropic cases. Expansion of the dimensionless bare coupling constants u on and the normal- 
ization constants ^( n ), Z^^ = Z^^Z^^) as functions of the dimensionless renormalized 
couplings u n up to two-loop order as 

u on = u n (l + a ln u n + a 2n u n ) , (3a) 
Zftn) = 1 + b 2n u 2 n + b 3n ul, (3b) 

Z^{n) = 1 + Ci n U n + C 2n u\, (3c) 

along with dimensional regularization will be sufficient to find out all critical exponents. 



III. SCALING THEORY FOR THE ANISOTROPIC CASES 

The anisotropic behaviors are characterized by correlation lengths £i, When con- 

sidered independently they define independent renormalization group transformations along 
the several competing directions. In momentum space, they induce independent flows in 
each external momenta scale k±, kl. 

In order to define the renormalized vertex parts we consider a set of cutoffs Aj (j = 
1, ...,L), each of them characterizing a different competing subspace. As functions of the 
bare vertices and normalization constants they read 

r £TW)> Qi(n), 9n, «») = 4n) Z ?(n)( r(JV ' L) (Pi(n), Qi(„), An, An) (4) 

— Sn,qSl,2^^ (Q(n),Q(n),K, A n ) Iq^ )=k 2 ) 

where p^ n ) (i = 1, N) are the external momenta associated to the vertex functions 
with N external legs and Qi( n ) (i = 1,...,M) are the external momenta associated to the 
M insertions of <fi 2 operators. From the last section, -u 0n , Z^i n \ and Z$ii n \ are represented 
as power series in u n . In order to write the renormalization group equations in terms of 
dimensionless bare and renormalized coupling constants, we shall discuss the central idea 
which underlies the subsequent scaling theory. 

Consider the volume element in momentum space for calculating an arbitrary Feynman 
integral. It is given by d^^^^qll^d^k^y Recall that q represents a (d — X^2 m *)" 
dimensional vector perpendicular to the competing axes and ku\ denotes an m^-dimensional 
vector along the ith competing subspace, respectively. The Lifshitz condition 5 0n = r nn r = 
suppresses the quadratic part of the momentum along the m 2 competition axes, the quadratic 
and quartic part of the momentum along the m 3 competing directions, and so on, such that 
the mi competing subspace is represented by a 2Lth power of momentum in the inverse 
free critical (t = 0) propagator, i.e., 1 {q,k) = q 2 + Xln=2 a n(tiL\) n - In order to be 

dimensionally consistent, the canonical dimension in mass units of the various terms in the 
propagator should be equal. 

Our normalization conditions give us a hint that we can get rid of the a n parameters pro- 
vided we make simultaneously dimensional redefinitions of the momenta components along 
each type of competition subspace in a complete analogy to the second character case. Let 
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[q\ = M be the mass dimension of the quadratic momenta. Since all momentum terms in 
the propagator should have the same canonical dimension, this requires that [ku\] = M<. 
These simultaneous dimensional redefinitions of the momenta along the competing axes 
is only possible due to the Lifshitz condition. The volume element in momentum space 

d a ~^i=2 mi qllf =2 d mi k^ has mass dimension M 1 . The dimension of the field can 

be found from the requirement that the volume integral of the Lagrangian density (1) is 

dimensionless in mass units. In other words, one obtains [(f)] = M2 1 L -= 2 * ' . In mo- 
mentum space the one particle irreducible (1PI) vertex functions have canonical dimension 



[T^(k l )]=M N+ ( d -^- 

Let us describe the theory in terms of dimensionless parameters. As the coupling con- 
stants are associated to T^ 4 \ we can write g n = u n (K%*)~2~ , and X n = Uo n (^ n )" 2 ~, where 
€l = 4 + Xln=2 ( ' n ~ 1 ^ m n — d. Expressed in terms of these dimensionless coupling constants, 
the renormalization group equation can be cast in the form: 

The functions 

Pn = ( K ™|)r)> ( 6a ) 

l4>{n)\U n ) = Pn 7^ (6b) 



n 



7^(n)(«n) = ~Pn ^ (6c) 

are calculated at fixed bare coupling X n . The /? n -functions can be rewitten in terms of 
dimensionless quantities as 

dlnu 0n , 1 

Pn = -ne L (— ) . (7) 

Note that the beta function corresponding to the flow in n n has a factor of n compared to 
that associated to the flow in K\. 

For the anisotropic case, the multi-parameters group of invariance is manifest in the 
solution of the renormalization group equation, which is given by 

r Rl)( k i(n),Un,Kn)=exp[~ f 7^( n ) (u n (p n )) — ] F^K , U n (p n ) , K n p n ) . (8) 

^ Jl x n 

From the above analysis, the dimensional redefinitions of the momenta along the distinct 
competing axes result in an effective space dimension for the anisotropic case, namely, (d — 
^2n=2 <n ~ 1 " > m n) • We discover the following behavior for the 1PI vertex parts under 
flows in the external momenta: 



, ., .'. ., ,■ \ - „'"" ' v " - 

N r Pn , , ^dx r 



r (JV) (o ki \ u k)- n ^ d -^^ m ^- 2 " — ) (Q) 



exp[-— / l<t,(n){u n (x n ))— ] 

* Jl x n 



U n(Pn), K n ) . 
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The behavior of the vertex functions at the infrared regime is worthwhile, since their fixed 
point structure will determine the scaling laws and the critical exponents for arbitrary m n - 
dimensional competing subspace. These L independent fixed points are defined by f3 n {u n ) = 
0. At the fixed points the simple scaling property holds 

1 R(n)v" n/ M n )' U "i K n) — Pn {LV) 

For iV = 2, we have 

^R(n)(Pn^(n)' U ni K n) — Pn H ^ n ^T^ R ^(k( n ), U* n , K n ) . (11) 

The quantity 7</,( n ) (?•**) can be identified as the anomalous dimension of the competing sub- 
space under consideration. 

This can be readily generalized to include L insertions of 2 operators such that the RG 
equations at the fixed point lead to the solution ((N, M) 7^ (0, 2)) : 

n\N+(d-^ L ("-^ m ) "^-^=2^^™™) 2tfl ^(n) , , f » 

Y {N ' M) (ok-< ^ an; n 77* k \ -n L "= 2 " } 2 1 2 +M Vw 

1 R(n) \P K i(n), PPi(n), «„, «nj — Pn 

X ^ii(ra) \^i{n) i Pi{n) ) K n ). (12) 

Writing this at the fixed point as 

r {N,M)( t, \_ "K'i-EL2 1 V 1 " ! »)-' v ' i *]+ M ^ (iV,M),, * \ n o\ 

1 R(n) \P K i(n), PPi(n),U n , Kn) — Pn 1 R(n) \ K i{n) , Pi(n) ; \i J , U<Jj 

the anomalous dimensions of the insertions of </> 2 operators are 0^2 = —In + 7^2^ (tt*). 

The scaling relations can be found by going away from the Lifshitz critical temperature 
(t 7^ 0) staying, however, at the critical region 5o n = r nn r = 0, which is the generalization 
of that from the ordinary second character Lifshitz behavior. Above the Lifshitz critical 
temperature, the renormalized vertices for t 7^ can be expressed as a power series in t 
around the renormalized vertex parts at t — 0, as long as iV 7^ 0. Then one can show that 
the RGE for the vertex parts when t 7^ are given by 

d d 1 d 

The key property of the solution is that it is a homogeneous function of the product of ku n ) 
(to some power) and t only at the fixed point 77*. As the value of u n is fixed at u n , we shall 
omit it from the notation of this section henceforth. Thus, at the fixed point the solution of 
the RGE reads 



r^ ) (fci(„),t,«n) = «n 2 F^\k i{n) ,K n t^). (15) 

Defining 6 n = —^2^, and using dimensional analysis, it is easy to show that 
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_R(n)V > «/ ~~ P n n 

F^iPn'hin^iPn 1 ^- 2 ^). (16) 

The choice p n = ft„(^) e ™+ 2n , can be substituted back in (16), implying that the vertex 
function depends only on the combination k^ n ^ n apart from a power of t. Since the cor- 
relation lengths £ n are proportional to t~ Un , it implies that the critical exponents v n satisfy 
the identity 

u- 1 = 2n + e: = 2n- 1 ; Hn) . (17) 
For convenience we could have defined the function 

7^(n){Un) = ~Pn 7^ • (18) 

In that case we would have discovered the equivalent relations 

v~ x =2n-7] n - 7 02(n) «). (19) 

For iV = 2 we choose p n — ^(n)-, the external momenta. Then ^\J n t, Kn) — 
k 2n ~ r,n n r ^- f{k(n)£, n )- The infrared regime corresponds to £ n — > oo and fy n ) — > such that 
f{k{ n )£,n) -> Constant. The definition /„ = {k^n) 2n ~ Vn f{k( n )^n), leads to 

r flW^(n),*,«n) = (^ ) a) 2n -' ? "C/n(A;(„)a). (20) 

(2) —1 

Since the susceptibility is proportional to t~ 7n as fc(„) — > 0, and T R ^ = X( n )i the suscepti- 
bility critical exponents are given by 

-f n = v n (2n-r] n ). (21) 

We now discuss the scaling law appropriate to relate the specific heat critical exponent 
to the others critical indices. The analysis of the RG equation for above T L at the 

fixed point yields information about the specific heat exponents. In that case it reads 

+ 7fc (B) (2 + *|))r§3 = K 2n )^B n (u* n ), (22) 

where B n (u^) is given by 

K 2n )^B n (u:) = -^W^^wWw;-^).^)!^- (23) 



The general discussion given up to now for the vertex part will be useful to uncover 

the homogeneous part of the solution. In fact, at the fixed point the generalization of the 
solution for is written as 

R(n) 



^S\Hn),Qi(n),t, K n ) = Kl U( " r * 2(n) Fi N > M \p i{n) ,Q i{n) ,K n t^). (24) 
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At the fixed point, the temperature dependent homogeneous part for h has the 
following property 

-2 * -i 
r ^),h(Q(n),-Q(n),t,K n ) = K n ^ W (Q (n) , -Q (n) , kJ*^) . (25) 

This is going to be identified with the specific heat at zero external momentum insertion 
Q(n) = 0. Using the dimensional analysis results, one can show that 

^ hh (Q(n),-Q(n ) ,t^n)=Pn ^ " V( " (26) 

X r S,h(Pn 1( 5(n), -^ 1 Q(n),P^ 2n *,Pn 1 «n) > 

and substituting this into the solution at the fixed point, it yields 

^R(n),h(Q( n )' ~ Q{n), t, K n ) = pn * ( ) Kn * ( ' (27) 

-l 

x F^ 2 \p^Q {nh -p- l Q {n) ,p- l K r Xp- 2n ty^)). 

Once more, choose p n = K n (-^) <?n+2n . Replace this in last equation, take the limit Q(n) - ► 
and identify the power of t with the specific heat exponent a n . The result is 

a n = 2 - n(d - V ^ - V )y„. (28) 

n=2 

The inhomogeneous part can now be discussed. Take — an d choose a particular 
solution of the form: 

C p (u n ) = (C)^(7 P K). (29) 
When this is replaced into the RG equation for at the fixed point, we learn that 



C p (u* n ) = (O * r < u B n (<)- (30) 

U n IL\U Z^n=2 n " Ln > Z 

Summing up both terms gives the following general solution at the fixed point: 

rg$ = K- 2 ")^^)- " + Vd £»«))■ (3i) 

The situation for T < Tl is as follows. For simplicity consider the case of magnetic 
systems. The renormalized equation of state furnishes a relation between the renormal- 
ized magnetic field and the renormalized vertex parts for t < via a power series in the 
magnetization M, i.e., 

H {n) (t,M,u n ,K n ) = w MNT{ it) } (hn) = 0;*,«n,«n), (32) 

JV=1 
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where the zero momentum limit must be taken after performing the summation. The mag- 
netic field satisfies the following RG equation: 



d „ d 1 T „ , d . d 



[Kr. 



+ &7i H N lm^ N + M 7w) + l^(n)tT-)H {n) (t, M, u n , n n ) = 0. (33) 



n dn n ' , n du n 2 mn)y ' dM J ' "rw"dt' 
The equation of state has the following form at the fixed point: 



H (n) (t, M, K n ) = K n 2 h ln (K n M^,K n f* 2 (") ). (34) 

Dimensional analysis arguments can be used to determine how a flow in the external mo- 
menta affects the renormalized magnetic field. The flow produces the following expression: 

H {n) (t, M, n n ) = p n n [ - + — L (n- 1)m) . -f)- (35) 

Pn 2 

The standard choice corresponds to p n being a power of M 

Pn = U r^r^ ] n[{d -^ ^^)-2 ]+ , n (36) 

rvn 

Replacing this into (35) and from the scaling form of the equation of state H^(t,M) = 
M Sn f(^^), we obtain the remaining scaling laws 

MS 

Pn = 7}i'n(n(d-^2 — . — nn) - 2n + i] n ) , (37a) 

1 i=2 ^ 



n(d ~ Ef=2 V 1 ^) + 2n - ?] n 

6n = (A ^ (i-D ^ 9 I ' (37b) 

n ( d ~ l^i=2 —r m i) ~2n + r] n 



which imply the Widom j n — (3 n (5 n — 1) and Rushbrooke a n + 2(3 n + j n = 2 relations for 
arbitrary competing or noncompeting subspaces, since n = 1,...,L. We note that there is 
one set of scaling relations for each competing subspace. This suggests that all the critical 
exponents take different values in distinct subspaces. We are going to see that this is not 
necessarily true, since the fixed point structure restricts the values of most critical exponents 
to be the same in different competing subspaces. It is a direct consequence that there is 
only one fixed point independent of the space directions under consideration. 

The perturbative calculation of the critical exponents and other universal quantities 
follows from a diagrammatic expansion whose basic objects are Feynman diagrams. We 
shall use the loop expansion for the anisotropic integrals with the perturbation parameter 
€l = 4 + X^n=2 ( '"~ 1 ' > m -n — d. The solution of the Feynman diagrams in terms of results 
in the e^-expansion for the universal critical ammounts of the anisotropic criticalities. The 
anisotropic integrals are described using the generalized orthogonal approximation in Ap- 
pendix A. This approximation yields a solution which is the most general one compatible 
with the homogeneity of the Feynman integrals for arbitrary external momenta scales. With 
this technique all the critical exponents in the anisotropic cases can be obtained as will be 
shown in Section V. 
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IV. SCALING THEORY FOR THE ISOTROPIC BEHAVIORS 



To begin with let us promote a slight change of notation with respect to the conventions 
presented in our previous letter [15]. There, the subscript associated to each type of m n 
competing axes was chosen as An. Here, we choose the subscript n to express the same 
thing. This will cause no confusion to the reader since the anisotropic and isotropic cases 
are considered separately in this work. Then an arbitrary ammount A 4n , should be changed 
to A n . In particular, the perturbative parameter discussed in section II is now represented 
as e n . Obviously, whenever d = m n , the volume element in momentum space is given 
by d mn k. Setting a n = 1 we perform a dimensional redefinition of the momenta such 
that [k] = M«. Accordingly, the volume element has dimension [<i mn /c] = . The 
dimension of the field in mass units is [<p] = M^" -1 . The 1PI vertex parts have dimensions 
[r^(/c„)] = M N+ ~^~ N ^ . Then, make the continuation m n = An — e n . The coupling 
constant has dimension A4 n = M sr^ = M^r. In terms of dimensionless quantities, one has 
the renormalized g n = u n (n n n )^ and bare A n = UQ n {K n n )^ coupling constants, respectively. 
Again, the functions 

fin = (^f^r 1 ) ( 38a ) 

dlnZftn) 

7*(n)(«n) = Pn—r^ (38b) 



dlnZ^ (n) 

! { n){U n ) = -fi n ^ (38c) 



are computed at fixed bare coupling constant X n . The beta functions in terms of dimen- 
sionless quantities are given by (3 n = —e n ( 9l Q^ 0n )~ 1 . Notice that the beta function for the 
isotropic case does not possess the overall factor of n present in the anisotropic beta function 
fi n obtained from renormalization group transformations over the m n -dimensional compet- 
ing subspace. This is a very close analogy to the second character behaviors and a general 
property of Lifshitz critical behaviors. 

The dimensional redefinition of the momenta along the m n competing axes leads to an 
effective space dimension for the isotropic case, i.e., (— ). Under a flow in the external 

momenta we find the following behavior for the 1PI vertex parts T^^: 

rfUPnk, u n , K n ) = p n n [N+ ^- N ^ ] expi-% JT i m {u n {x n ))^\ (39) 

X ^i?(n)(^' U n(pn)j K n) ■ 

Note that since there is only one type of space directions in the isotropic behaviors, we 
do not need to use a label in the external momenta specifying the type of competing axes 
considered as we did in the anisotropic cases. At the fixed point, the simple scaling property 
for the vertex parts follows: 

C(M„ <■ O = ^+*-"*'-?WO (40) 
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For N = 2, we have 

^R(n)(Pnk> U ni K n) = Pn ^ ' ^fl^ (fc, W* , « n ) . (41) 

In the noninteracting theory d^ = -° 1 is the naive dimension of the field. At the 

m 

isotropic fixed point, the presence of interactions modify it such that — -f - — 1 + ^. The 
generalization to include L insertions of 4> 2 operators can be written at the fixed point as 
((iV,L)^(0,2)): 



T {N ' L) (o k- o v- v* k ) - o " 2 2 L V(n)p(iV,L) r , * v (A2) 

R(n) yPn^ti PnPi, U n , tv n ) — fJ n L R(n) "niM' 



Writing at the fixed point 

T {N ' L) (o k o n v* k ) ~ ™"-^ +i V r Wi)/L * n / 4 ox 

the anomalous dimension of the insertions of ^ 2 operators is given by e^2 = — 2n + 7^2 („)(«*). 
Above the Lifshitz critical temperature we find the following RGE 

1 

][Kn 'dK~ + l3n du' ~ 2 N7<i>{n)<KUn ^ + ^H^iu^t^T^ik^t,^, K n ) = °- ( 44 ) 

The solution at the fixed point is given by 

r?^'W.^>n) = k^F^,^^). (45) 
If we define 9 n = — 7^2 ( n ), we can use dimensional analysis to obtain 

X ^S ) (Pn 1 ^.(P» 1 «»)(Pn 4 *)^)- (46) 

We can choose p n = K n {-^:) 6n+2n , and replacing it in the last two equations, the vertex 
parts depend only on the combination ki£ n apart from a power of t. As £ ra is proportional 
to t~ Vn we can identify the critical exponent v n as 

u- 1 = 2n + e* n = 2n- 1 ; 2{n) . (47) 

For the sake of convenience we define the function 

_ , x n dln{Z^ {n) Z m ) 

l^(n){U n ) = -Pn 7^ • (48) 

In terms of this function the last equation turns into the following relation 

V n 1 = 2n - Vn ~ 7^(n)«)- (49) 

For iV = 2we choose p n = k, the external momenta. When £ n — > 00 and — * 0, simulta- 
neously, then f(k£ n ) — > Constant. The susceptibility is proportional to t~ 7n as fcj — > 0. As 
= the susceptibility critical exponent follows 
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7n = ^n(2n - T) n ). (50) 



From the RG equation for above T L at the fixed point, the scaling relation for the 
specific heat exponent can be found. The RG equation is 



where 



:< 2n )-5n«) = -^(n)^^-[r^ 2) (Q;-g,A n )| Q2 = K d- (52) 



The can be generalized to 



1 Nil^-L-r: 



fixed point as 



The homogeneous part of the solution for h is temperature dependent and scales at the 



-27 



rg^CQ.-Q,*,^) = «n * (n) Fi°' 2 \Q,-Q,K n f^). (54) 

This vertex function is going to be identified with the specific heat at zero external momen- 
tum insertion Q — 0. Use of dimensional analysis yields the result: 

rg$, fc (Q, -0, «n) = ^^^r^fa'Q, -p?Q, Pn 2r % Pn^n)- (55) 

Substituting this equation in the solution at the fixed point leads to 

n rmn_4]+2 7 * -2 7 * 



(56) 

The choice p n = K n {-^) e "-+ 2n can be made. Substitution of this choice into the last equation 
in the limit Q — > and identifying the power of t with the specific heat exponent a n , we 
obtain: 

a n = 2 - m n f„. (57) 

The inhomogeneous part can be found by taking Q = and choosing a particular solution 
in the standard way. Therefore, the general solution at the fixed point is given by 

Next, let us concentrate ourselves in the scaling relations when the system is below the 
Lifshitz critical temperature T < T^. The relation among the renormalized magnetic field, 
the renormalized vertex parts for t < and the magnetization M is given by 
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OO j 

tf (n) (t, M, u n , n n ) = t^M^JJ^^ = 0; t, u n , (59) 



TV! 

JV=1 



The RG equation satisfied by the magnetic field is: 

The solution of the equation of state at the fixed point has the following property 
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H {n) {t,M,K n ) = K n 2 h n {K n M^,K n f^)). (61) 

The scale change in the magnetic field followed by a flow in the external momenta can be 
written in the form: 



rill". 



H {n) (t, M, n n ) = p^ +1] H n (-L, ^tt v -)• (62) 

Pn p n l2n 1 Pn 

The flow parameter p n is chosen as a power of M such that: 

Pn = k„ [ [2 ff_ n] ] ^-L+ Vn j (63) 

and from the scaling form of the equation of state H( n -)(t, M) = M 5n /( n )(— t —), we obtain 
the following scaling relations: 

m n + 2n-r] n 

$n = - , (64a) 

m n -2n + r] n 

Pn = 2 U n(m n - 2n + r) n ), (64b) 

which imply the Widom j n = /3 n (5 n — 1) and Rushbrooke a n + 2/3 n + 7„ = 2 relations. 

Except for some minor modifications, the renormalization group treatment of each 
isotropic behavior is equivalent to treat separately each competing subspace appearing in 
the most general anisotropic behavior. It is easy to see that all scaling laws reduce to those 
from the usual critical behavior described by a (ft 4 field theory for n — 1. Then, the usual 
critical behavior is actually a first character isotropic Lifshitz critical behavior. For n — 2, 
they easily reproduce those associated to the second character Lifshitz point [8]. Therefore, 
a new nomenclature emerges from the study of these higher character Lifshitz critical be- 
haviors: the number of neighbors coupled through competing interactions is a fundamental 
parameter, generalizing the concept of universality class. Thus, the universality classes of 
isotropic behaviors (d = m n ) are characterized by (N,d,n). These statements will be put 
on a firmer ground when we calculate the critical exponents, as we shall see in the next 
sections. 

The isotropic behaviors are calculated using the generalized orthogonal approximation as 
well as exactly, without the resource of any approximation. The approximation is useful to 
complete the unified analytical description of the higher character Lifshitz critical behavior in 
its full generality, at least at the loop order considered here. On the other hand, the analytic 
exact (perturbative) solution is a conceptual step forward towards a better comprehension of 
this sort of system. At this point, the reader should consult the Appendix B in order to access 
the computation of the Feynman integrals using the generalized orthogonal approximation 
and Appendix C to see the exact computation required to finding the critical exponents. 
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V. CRITICAL EXPONENTS FOR THE ANISOTROPIC BEHAVIORS 



In this section we compute the critical exponents using the generalized orthogonal ap- 
proximation with the results derived in Appendix A. First, we attack the problem using 
normalization conditions. Second, the results are checked using a variant of the minimal 
subtraction scheme first developed in [8] for the anisotropic second character m 2 -fold Lifshitz 
point. 



A. Normalization conditions and critical exponents 

The bare coupling constants and renormalization functions were defined in section II. 
They are given by 

u on = u n (l + a ln u n + a 2n u 2 n ), (65a) 

Zct>( n ) = 1 + b 2n u 2 n + b 3n iil, (65b) 

^0 2 (n) = 1 + c ln u n + c 2n ul, (65c) 

where the constants a,i n , bi n , Ci n depend on Feynman diagrams computed at suitable symme- 
try points. The critical exponents associated to correlations perpendicular or parallel to the 
arbitrary competing m„-dimensional subspace can be calculated through the specification 
of the corresponding symmetry point. 

In terms of the constants defined above, the beta functions and renormalization constants 
can be cast in the form: 

Pn = -ne L u n [l - a ln u n + 2{a\ n - a 2n )u 2 n ], (66a) 
7<«n) = -ne L u n [2b 2n u n + (36 3 „ - 2b 2n a ln )u 2 n ], (66b) 
7 02(n) = ne L u n [c ln + (2c 2 „ - c\ n - a ln c ln )u n ]. (66c) 

The above coefficients can be figured out as functions of the integrals calculated at the 
symmetry points. We find 

ain = N „ + 8 [1 + h mL e L ], (67a) 
oex, 

{ N + 8 2 (N + $l (3AM-14) 1 

b 2n = -^Jf t 1 + ( 2 >W + \)e L ], (67c) 
b - ( N + L2 )( N + 8 ) | (N + 2)(N + 8) I 1 

hn ~ 12964 IOSTl [ ~l mL + 48 } ' (6?d) 

(N + 2) 

c\ n = — r [l + h mL e L ], (67e) 

6e L 

(iV + 2)(7V + 5) (jy+^Cg + j) i 

36e| 3e L 1 3 L 4 J K J 

These expressions are sufficient to find out the fixed points at 0(e|), which are defined by 
Pn( u n) = 0- As was seen in the last section, every integral computed at arbitrary symmetry 
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points SP\, SPl gives the same result, irrespective of the considered subspace. The 
overall factor of n = 1, ...,L in the /3 n functions drops out at the fixed points, such that 
the renormalization group transformations realized over and kl will flow to the same 

fixed point given by (u\ — ... — u* L = u*) 



(9N + A2) 
+ ( 8 + iV)2 



(68) 



It is instructive to separate explicitly the noncompeting and competing subspaces. The 
functions 7^(1) and 7</>2(i) read 

(iV + 2) r , , L n 9 (N + 2)(N + 8) o , « 

7«i) = ^2 t 1 + ( 2/ W + 4)^]^ - ^ ^ i' (69) 

7</> 2 (i) = ^— ^— ^it 1 + h mL e L - ifzi]. (70) 

When the value of the fixed point is substituted into these equations, using the relation 
among these functions and the critical exponents f]i(= r] L2 ) and V\(= Vl2)-> we find: 

1 2 N + 2 M ,6(3iV + 14) 1,. 

* = 2 €l WTW [ + £L( (iV + 8^ " 4 )] ' (71) 

1 (JV + 2) l (iV + 2)(iV 2 + 23iV + 60) 2 

2/1 " 2 + 4(ivT8) eL + 8 (N + W ' (72) 

The coefficient of each power of 6l is the same as that coming from the second charac- 
ter Lifshitz behavior m 3 = ... = mi = 0. Consequently, the reduction to the Ising-like 
universality class 777-2 = case is warranted. The several beta functions corresponding to 
distinct competing axes satisfy the property j3 n = n(3\. This implies that J^n) = n l^>(i) an d 
7</> 2 (n) = ^702(1) • Then, we have 

n (N + 2) 6(37V+14) 1 

^ = 2 (ivT^F [1 + ei( (N + sy - (73) 

11 (yY + 2) l (iV + 2)(iV 2 + 23iV + 60) 2 
^ = n { 2 + W^) L 8 (JVT8p (74) 

At 0(e 3 L ), the relation 77 ra = nr\\ is satisfied, whereas at 0(e 2 L ), the relation z/ n = -v\ holds. 
Strong anisotropic scale invariance [24] is exact to the perturbative order considered here, 
and within the generalized orthogonal approximation it is expected to hold at arbitrary 
higher loop order. Differently from the critical indices rj n and v n which depend explicitly on 
the m^-dimensional subspace under consideration, the other exponents take the same value 
in each subspace even though they are obtained through independent scaling relations along 
the distinct competing axes. They are given by 

(N + 2) (N + 2)(N* + 22N + 52) ^ 2 

il = 1+ w^) L — w^w — (75) 

_ (4-AQ (N + 2)(N* + 30N + 56) 2 

aL - w+8) L w+w — (76) 
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2 2(iV + 8) 



1 3 



e L + 



(7V + 2)(2iV+l) 
2(iV + 8) 3 



(77) 



h = 3 + e L + 



(N 2 + UN + 60) 
2(iV + 8) 2 



(78) 



The exponents correctly reduce to those from the second character behavior [8,15], with 
a further reduction to the Ising-like case when m 2 = 0. In fact the universality classes 
reduction of generic higher character anisotropic Lifshitz points to that from Ising-like critical 
points is manifest in all critical exponents. Hence, this universality class reduction is a 
generic property of arbitrary competing systems. 

To check the correctness of these exponents, it is convenient to calculate them in an- 
other renormalization procedure. Let us check these results using minimal subtraction of 
dimensional poles, as we are going to show next. 



In the minimal subtraction renormalization scheme, the common situation is to have 
just one momenta scale /x [25] which is called k in the present paper. Nevertheless, the di- 
mensional redefinitions performed over the external momenta characterizing arbitrary types 
of competing axes permit a picture of the anisotropic cases with L independent momenta 
scales. 

The calculation of the critical exponents along an arbitrary kind of competition subspace 
can be done, provided all external momenta not belonging to that subspace are set to zero. 
Then, we define Kj to be the typical scale parameter of the jth subspace, calculate the 
renormalization functions for arbitrary external momenta along the m^th space directions 
and require minimal subtraction of dimensional poles. This procedure is inspired in the 
method formerly discussed in the second character anisotropic Lifshitz behaviors [8]. Thus, 
although the external momentum associated to the competing subspace under consideration 
is kept arbitrary in all stages of the computation, the same is not true for all other external 
momenta corresponding to distinct competing subspaces. This restriction on the values of 
all the external momenta is the price to be paid in order to describe independently the 
scale transformations of each inequivalent subspace. This is the main difference of this min- 
imal subtraction scheme with several independent momentum scales from the conventional 
method with just one momentum scale. 

Here we will content ourselves in showing that the diagrammatic procedure to calculate 
the fixed point using minimal subtraction results in the same functions j^n) and 7^2 ( n ) at 
the fixed point as those coming from normalization conditions. This is equivalent to prove 
the renormalization scheme independence of all critical indices. 

In minimal subtraction, the dimensionless bare coupling constants and the renormaliza- 
tion functions are defined by 



B. Minimal subtraction and critical exponents 




(79a) 



i=i 



oc 




(79b) 



i=i 
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^02(n) = 1 + C in{tL)u % n . (79c) 
i=l 

The renormalized vertex parts 

^i?(n)(^(™)' W ™' = ^</>(n)r| n ^(A; n , Won; K n); (80a) 

^ii(n)(^(")' ' U ™' K ") = ^'|(n)'' n (n)(^ i (™)' W 0™' (80b) 
r R$)( k l(n),k 2 (n),P(n);U n ,K n ) = Z^^T^ (fc 1(n) , fc 2 ( n ), p (n ) , M „, K n ) , (80c) 

are finite by construction when — > 0, order by order in u n . One should bear in mind 
that the external momenta in the bare vertices are mutiplied by k" 1 . Since = pi are 
the external momenta perpendicular to the competing axes, whereas k^ n ) = k' i( ^ are the 
external momenta parallel to the m n -dimensional type of competing subspace, the coeffi- 
cients aj n (ei), bi n (ei) and Ci n {e£) are obtained by requiring that the poles in be minimally 
subtracted. The bare vertices are written in the form 

r£)(fc(n), U 0n, «n) = fc^(l - B 2n ul n + B 3n U 3 0n ), (81a) 
r (n)( fc i(")' W 0n, «(n)) = <>0n[l ~ ^ln^On + (^2n + -42) M 0rJ> ( 81b ) 
^'^(^(n), k 2 (n),P(n); U 0n , K n ) = 1 - C ln U 0n + (C^ + C^wL ( 81c ) 

Remember that B 2n is proportional to the integral /3 and B 3n is proportional to I 5 which 
are calculated with all external momenta not belonging to the m n -dimensional subspace set 
to zero. 

The coefficients are expressed explicitly by the following integrals: 

A ln = ( N + 8 ) [l2{ k Hn) + k 2{n)) + h{ k m + k m) + h{ k 2{n) + k 3{n))i (g2a) 

18 K n K n K n 

A W _ + 6iV + 20) 2 fc 1(TO) + fc 2(TO ) 2 ^ 1( n)+fc 3 (n) ^ 7 2( fc 2(n)+fc 3 (n) ^ ( ^ 

108 K n K n K n 

4 2 J = * [i4(— ) + 5 permutations], (82c) 
54 K ra 

B 2 „ = (82d) 
18 K n 

= (AT + 2)^ + 8) fa 

3n 108 y n n h v ; 

C7 ln = ^±l [h{ k nn) + k 2{n) ) + j2( fc M± A^ ) + /2( ^) + ^ ) )L (g2f) 

= (jV + 2) 2 2 jlH + fc 2(TO ) 2 fc 1(B) + fc S(B ) 2 jg(n)+j3(n) ( } 

108 /t ra /t n /t n 

C£ } = ^_t^[j 4 (^M) + 5 permutations}. (82h) 
36 re n 

We have at hand all we need to determine the normalization constants at least at two-loop 
order. Requiring minimal subtraction for the renormalized vertex parts above listed, it can 
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be verified that all the logarithmic integrals depending upon each arbitrary (nonvanishing) 
external momenta subspace appearing in / 2 ,-^3,^4, and 1$ cancell out. This leads to the 
following expressions for the normalization functions and coupling constants: 



UQt. 



M (AT + 8) r (iV + 8) 2 (3N + U). 2 , 
u n (l + ^ -u n + [ v — J - 'Aul 



J <f>(n) 



1 - 



= 1 + 



6e L 
N + 2 2 

T^ Un + [ ~ 

N 



36e| 



24e, 



(iV + 2)(iV + 8) (N + 2)(N + 



1296e| 



+ 



5184e L 



\u r 



6e ; 



^ + [ (iV+ S ( f +5) -^K). 



36e| 



24e, 



Using the renormalization constants we obtain: 



AN + 2) 2 (iV + 2)(7V + 8) 31 



l<t> 2 {n) = 



71- 



72 

(N + 2) 

6 



1728 



1 . 
u n [l - -u n \ 



(83a) 
(83b) 
(83c) 

(84) 
(85) 



The fixed points are defined by /3 n (u*) = 0. Then, we learn that the fixed points 
generated by renormalization group transformations over Ki,..., and kl are the same, namely 



6 



u. 



% + N 



e L < 1 + e L 



(9A^ + 42) 
(8 + N) 2 



(86) 



When this result is replaced in the renormalization constants at the fixed point it yields 
7^(„) = Vn, where rj n are given by Eqs. (71) and (73). In addition, we have 



7>M = U 



{N + 2) 
(N + S) 



e L [l + 



6(N + 3) 
{N + 8} 



(87) 



The reader can verify that the exponents v n encountered by using the last equation are the 
same as those obtained via normalization conditions Eqs. (72) and (74). This proves the 
consistency of either renormalization scheme for the anisotropic Lifshitz critical behaviors. 



C. Discussion 



The generic Lth character Lifshitz critical behavior naturally extends the comprehen- 
sion of the usual second character Lifshitz criticality. The latter is characterized by one 
noncompeting subspace and only one competing subspace, whereas the former is character- 
ized by several types of competing subspaces. The expressions for the critical exponents can 
be analysed to extract further information concerning those systems. For instance, when 
m,3 7^ and 777,2 = m4 — ■■■ — % = 0, the third character behavior is recovered and correctly 
reduces to the Ising-like behavior for 7773 = 0. The main characteristic of generic third char- 
acter behavior is that there are m 2 ,m 3 7^ competing axes with m 4 = ... = tul = 0, and 
so on. From a phenomenological perspective in magnetic systems, may be it is worthy to 
assemble all magnetic materials presenting Lishitz critical behavior and analyse their critical 
exponents. Choosing those alloys in the same conjectured universality class, the greater the 
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difference in their critical exponents the more likely they are in an alternative universality 
class contained in the CECI model analysed here. 

It is interesting to note that exact strong anisotropic scale invariance [24] is valid in 
the CECI model as a result of the generalized orthogonal approximation. Since the model 
describes the physics of short ranged competing systems, the limit L — > oo in the Lth 
charater Lifshitz point should not be taken, since it would describe long range competing 
systems. Nevertheless, since no restriction on L was made in the beginning of the discussion, 
the CECI model can be viewed as describing a particular type of long range competing 
interactions. If we go on and take this limit in the expression of the critical exponents we 
find that r] L tends to infinity, whereas — ► as a consequence of the strong anisotropic 
scale invariance. This implies that close to the Lifshitz critical temperature the correlation 
length £l does not diverge in that limit, instead of having a usual power law divergence. 
This fact is a nonperturbative result valid to all orders in perturbation theory within the 
context of the e^-expansion. 

Another feature emerges from this limit by looking at the critical dimension d c = 4 + moo. 
(Recall that in the Lth anisotropic character critical behavior the system has mi competing 
axes and d — rriL noncompeting space directions, i.e., rri2 = ... = tul-i = 0.) This limit yields 
a mechanism which is a natural new way to study extra dimensions without destroying the 
renormalizability of the corresponding field theory, while retaining the nontrivial aspect of 
the fixed point. In spite of the divergence of the anomalous dimension and the vanishing 
of the correlation length exponent along the ttil competing directions in this limit, all 
other exponents are well behaved and have the same value of those corresponding to space 
directions without competition, as can be seen explicitly by the expressions for the exponents. 
This is so since the anisotropic scaling laws only contain the safe combination Lvl which is 
always finite in the limit L — > oo. Once more, note that — > reduces to the previous 
(p 4 universality classes. Recall that the isotropic case is characterized by d — mi = 4L. 
Some care care must be taken in the interpretation of this case in the limit L — > oo, for 
the approach from the anisotropic (0 < tul < 4L) to the isotropic case takes infinite steps, 
which is not as simple as in the case when L is finite. This point deserves its own analysis 
for future work. 

On the other hand, the last few years have witnessed some ideas in quantum field theo- 
ries that resemble very much issues contained in the analysis of Lifshitz critical phenomena. 
The closest analogy with the Lifshitz field theoretic tools presented here is the very recent 
idea of ghost condensation producing a consistent modification of gravity in the infrared 
(long distance limit) [19,26,27]. In this framework, gravity is modified to have attractive as 
well as repulsive components. The latter mimics a kind of dark energy [28], whose ghost 
condensate is a physical fluid arising from a theory where a real scalar field changes with 
a constant velocity. The ghost condensate appears as the physical scalar excitation around 
the background generated by the scalar field whose vacuum expectation value is defined by 
a constant value of its time derivative. Consequently, the effective action for the field repre- 
senting the ghost condensate has kinetic terms with quadratic time derivatives and quartic 
space derivatives of the field, therefore breaking Lorentz invariance [19]. The instability in 
momentum space is manifest in the absence of kinetic terms quadratic in the space deriva- 
tives of the ghost condensate. This characteristic is a result of the competing nature between 
the attractive and repulsive components of the gravitational force. This is a precise analogy 
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with the Lifshitz critical region in the case of the usual second character behavior included in 
the discussion given in the present paper. The utilization of an analogous reasoning leads us 
to conclude that the CECI model can be used to extract further insights from these new ef- 
fective quantum field theories where Lorentz invariance is broken. It permits generalizations 
for the ghost condensate when higher powers in space derivatives of this field are present in 
its corresponding effective action. For example, at large distances the case of quadratic time 
derivatives and 6th space derivatives in the kinetic term would correspond to a gravitational 
interaction with attractive/repulsive/attractive competing components and so on. Further 
analysis of the model might be helpful to address the perturbative calculations regarding 
the ghost condensate in that scalar field background. 

Numerical methods to probe the results obtained here are in their infancy for the CECI 
model. Earlier Monte Carlo simulations for a uniaxial third character behavior were per- 
formed [11] and the existence of the corresponding Lifshitz point at nonzero temperature 
was established. Unfortunately, no critical exponent was determined for the third character 
Lifshitz point. Since then, perhaps due to the fact that these higher character criticalities 
were not well understood theoretically, these methods are still waiting for more investiga- 
tions. Recently, a model with antiferromagnetic couplings between nearest neighbors as 
well as antiferromagnetic exchange interactions between third neighbors was studied using 
Monte Carlo simulations and some quantum properties were investigated [29]. No ferro- 
magnetic couplings appear among arbitrary neighbors. Even though there is a quantum 
Lifshitz point there, it pertains to a universality class which might be different from that 
representing third character critical points discussed in the present paper. Further numerical 
studies motivated by the CECI model are necessary to understand completely the classical 
and quantum properties of its critical regions. 

VI. ISOTROPIC CRITICAL EXPONENTS IN THE ORTHOGONAL 

APPROXIMATION 

We now address the calculation of critical exponents using the generalized orthogonal 
approximation from the results obtained in Appendix B. As we have seen in the calculation of 
Feynman integrals, this problem can also be tackled without performing any approximation 
during all steps of the calculation. This technique shall be postponed until next section. In 
a way or another, the approach is equivalent to treat each competing subspace separately. 
Though very similar, the framework of this section turns out to be more economical than 
that reported in the anisotropic cases. 



The basic definitions of the bare coupling constants and renormalization functions were 
encountered before and are given by 



A. Critical exponents in normalization conditions 





Z^n) = 1 + b 2n u 2 n + b 3n Un 
Z(f,2(n) = 1 + c ln u n + c 2n u: 




24 



where the constants a,i n , bi n , C{ n depend on Feynman integrals at the symmetry point called 
SP n . Let K n denote the competing m n -dimensional subspace in the isotropic cases. 
The beta-function and renormalization constants can be expressed in the form 

Pn = -e n u n [l - a ln u n + 2{a\ n - a 2n )u% (89a) 

70(n) = -e„M„[26 2 nMn + (36 3 n - 1b 2n ain)u% (89b) 

70 2 (n) = £nUn[Cln + (2c 2n - c\ n - a ln C ln )u n }. (89c) 

The coefficients above obtained as functions of the integrals calculated at the symmetry 
point read 

din = N r + 8 [1 + 77- £ "L ( 90a ) 
be n In 

,N + 8. 9 r 2JV 2 + 23iV + 86 1 , 

tt2 » = ( ^r )+l — ^ — i- (90b) 

(iV + 2) r 5 . , , 

3 1296ne 2 5184n 2 e„ ' V ; 

(N + 2) r 1 , . . 

dn = — R -1 + ^-en, 90e 

6e„ 2n 



= (iV + 2)(^ + 5) (JV + 2)(2Af + 7) 
36e? 72ne, 



~n 



The equation /?«(«*) = defines the fixed point. Thus, we find 



6 L 1 



8 + AT n 



1 (9N + 42) 
~2 + (8 + iV) 



2 



(91) 



We stress that this fixed point is different from that arising in the anisotropic behavior and 
cannot be obtained from it within the e^-expansion described above. The functions 7^,( n ) 
and 7^2 ( n ) are found to be 

(iV + 2) M 5 . 2 (iV + 2)(iV + 8) 3 
(TV + 2) . 1 1 . 

7<^ 2 (n) = g ttn[l + 2^£n ~ 7^ U n\- ( 93 ) 

When the fixed point is replaced inside these equations, using the relation among these 
functions and the critical exponents r\ n and u n , we find: 

n = ^ N + 2 Fl + 6 1( 6(37V + 14) - 1)1 (94) 
^ n 2n 6 ™(iV + 8) 2[ + n n { (N + 8) 2 4 JJ ' 1 j 

J_ (iV + 2) i (AT + 2)(iV 2 + 23iV + 60) 2 

n "2n + 4n 2 (iV + 8) " + 8n3 (TV + 8) 3 n ' 1 J 

The coefficient of the e 2 term in the exponent r\ n is positive, consistent with its counterpart 
in the anisotropic cases as well as in the Ising-like case. In the generalized orthogonal 
approximation the competing momenta are not sufficient to induce its change of sign. 
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Now using the scaling relations derived for the isotropic case we obtain immediately 



(N + 2) (N + 2)(N^ + 22N + 52) ^ 2 

ln ~ + 2n(N + 8) n+ 4n 2 (iV + 8) 3 1 ' 

(4-iV) _ (iV + 2)(iV 2 + 3(W + 56) 2 
n 2n(7V + 8) n 4n 2 (iV + 8)3 

* 1 3 F i (iV + 2)(2iV + l) 2 

Pn 2 2n(iV + 8) n+ 2n 2 (iV + 8)3 

j -3 | \ | (^ 2 + 14iV + 60) 
dn " d + n Cn+ 2n 2 (iV + 8) 2 ^ 

The explicit dependence on the number of neighbors coupled through competing interactions 
is manifest in the above exponents. Hence, the universality classes for an arbitrary isotropic 
(d = m n ) competing system are determined by (N,d,n). The interesting fact is that the 
results for the ordinary critical behavior are correctly recovered in the limit n — > 1. 

The orthogonal approximation is a good approximation even for isotropic higher charac- 
ter Lifshitz critical behaviors for at least three main reasons. It preserves the homogenity of 
the Feynman integrals in the external momenta scales. Second, it indicates which param- 
eters are important to describe the universality classes of isotropic systems. And last, but 
not least, it manifests one of the most important properties of competing systems, namely, 
the reduction to the Ising-like universality classes in the limit when all interactions beyond 
first neighbors are turned off. 

It is important to emphasize a technical detail concerning the approximation just em- 
ployed. The leading singularities from the one- and two-loop diagrams contributing to the 
four-point 1PI vertex part do not get modified from the usual 4 . On the other hand, the 
leading singularities of the two- and three-loop for the two-point 1PI vertex function do get a 
factor of - with respect to those from the usual critical behavior. These integrals also do not 
change signs under the generalized orthogonal approximations. The calculation performed 
without approximations shows that the leading singularities of diagrams contributing to the 
1PI two-point function change in a more complicated way, also changing sign depending on 
the value of n. We shall compare the differences in the values of the exponents utilizing 
the orthogonal approximation and the exact treatment later on. In particular, we shall per- 
form a numerical analysis for the isotropic second character behavior to understanding the 
deviations in both approaches. 

In order to check these results, let us analyse the situation using the minimal subtraction 
scheme. 



B. Critical exponents in minimal subtraction 



Minimal subtraction of dimensional poles in the renormalized vertex T^L, can be used 
to show that all momentum-dependent logarithimic integrals are eliminated in the renor- 
malization process leading the bare dimensionless coupling constant to be written as 

(iV + 8) iiV + 8) 2 (37V + 14V 9l , . 
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The fixed point can be found out 



Un (N + 8) n + n(iV + 8)3 e ™' 



In addition, the normalization constants are given by: 

(iV + 2)^ 



(N + 2)(iV + 8) (N + 2)(N + 8) 

+ [ 1296ne 2 + 5184n 2 e n K ' 1 ° j 
- (TV + 2) 

A£ 2 (n) = 1 H u n 

Consequently, the functions ^(n) an d l<t> 2 (n) can be obtained directly 

(N + 2) 2 (iV + 2)(iV + 8) 3 . 

7<6fn) = - (104a) 

WW 72n ™ 1728n 2 n v ; 

7^(n) = (iV ^ 2) K - ^u\). (104b) 

Substitution of these expressions in the function at the fixed point, one obtains the 
value of rj n as obtained in (94). The function 7^ 2(n ) at the fixed point reads 

(N + 2) . 6(N + 3) . /1A . 

Vw = jw^f^ 1 + n(iV + 8 ) 2 e J' ( 105 ) 

that is the same as the one obtained in the fixed point using normalization conditions. 
Therefore, it yields the same critical exponent v n from (95) as can be easily cheched. This 
constitutes the equivalence between the two renormalization schemes. 



VII. ISOTROPIC CRITICAL EXPONENTS IN THE EXACT CALCULATION 

We now turn our attention to the calculation of the isotropic critical exponents exactly, 
i.e., not using the orthogonal approximation. The algorithm we need to employ to obtain 
the critical indices is pretty much the same as that used in the orthogonal approximation, 
since the renormalization program and the scaling laws are approximation independent. 
The difference is that one should replace the Feynman integrals by their exact values already 
calculated in Appendix C. As we are going to discuss the case n = 2 in the next section using 
normalization conditions and minimal subtraction, we shall take normalization conditions 
in this section. Nevertheless, with the resources furnished in the text, the reader should be 
able to check the exponents using minimal subtraction as well. 

Using the definitions given for the power series of the bare dimensionless coupling con- 
stant and renormalization functions in terms of the dimensionless renormalized coupling 
constant, we find the following values for the coefficients of each power of u n \ 
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a ln = ^-^[l + D(n)e n ], (106a) 
6e n 

N + 8 2 1 (AT 2 + 21AT + 86)£>(re) 57V + 22 (AT + 2)(-l)"r(2n) 2 
fl2n ~ ( 6e n ^ + 18 18 36r(n+ l)r(3n) 

.^(V'.^^),, (108b) 
p=i ^ p=i ^ 

, 1X „ (iV + 2)r(2n) 2 ri , 3 3^1 

hn = { - iy 72T(n + l)T(3n)e n [1 + ^ + 4 " 2 £ p 

1 n— ^ 1 3 ^ 1 

— + ^-) e J> ( 106c ) 

2 n — p 2 p 

p=l ^ p=3 ^ 

(iV + 2)(iV + 8)r(2n) 2 1 1 J(n) 1 
3n 1 J 108I> + l)r(3n) [ 6e 2 e„ l 3 24 

2rt— 1 n— 1 -. -. 3n— 1 

p=2 ^ p=l ^ p=3 y 

c ln = { -^-^-[l + D(n)e n ], (106e) 

" p=i p=i 

Here D(n) = |^(2n) —ip{n) + |"0(1). The fixed point at two-loop level is defined as the zero 
of the /^-function and it is given by the following expression: 

.... _ *. + ^, i ( (-i)-rw+4) + {20N + 88)(1 _ m) _ m 



n (N + 8) L L (A^ + 8) 2V r(n + l)r(3n) 

v ' p=l J p=l 

The renormalization functions 7</,( n )(wn) and 7</>2( n )(w n ) can be expressed in the following 
simple manner in terms of u n : 

, x / (iV + 2)r(2n) 2 . . . . 3 3 2 ^1 1 S=i 1 

7*.,(«.) = (-'r' ^+Drw ' 1 + (B(n) + 1 - 2 g ; + 1 g — p 

+ 2^p'* ' ' 216r(n+l)r(3n) '2 + ^p £f„-p ^ p K> (108) 

7^2(n)(«n) = -^y-^- Mn[l + ^(w)Cn ~ ^W«n]- (109) 

Substitution of the fixed point in the first equation gives directly the anomalous dimensions 
i] n . Using the scaling law relating the second expression with the exponents r\ n and u n , one 
obtains the exponent v n . Thus, we get 
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„ - { (^ + 2)r(2n) 2 (iV + 2)r(2n) 2 F(iV,n) 

/n 1 ; (A^ + 8) 2 r(n + l)r(3n) n 1 J (AT + 8) 2 I\n + l)T(3n) ' J 



where 



F(N, n) = [((_i)" r ( 2n ) 2 ( 4 ^ + 8 ) + (4QAT + 176)£>(ra)) ' 



r(n + l)r(3n) v 7 v 7/ (iV + 8) 2 

„ 2n— 1 ^ n— 1 ^ 3n— 1 ^ 

i-e-^iv^ 1 ' < m > 

p=i p=i p=i 



Analogously, 



n 2n 4n 2 (iV + 8) n + 4n 2 (iV + 8) 3 " U J 1 ' T(n + l)r(3ra) 2n 
+(14JV + 40)D(n)]. (112) 

Now, we use the scaling relations to obtain the remaining critical exponents. We find: 

. 1 i (iV + 2) , I (iV + 2) , 2 [r ir 2n(27V + 4)r(2n) 2 

7 " = 1 + 2n(iV + 8) e " + 4n 2 (iV + 8)3 e " [( - 1) T(n + l)r(3n) + + 2)(* + 8 ) 

+2w(14iV + 40)D(n)], (113) 

Qn " 2n(iV + 8) e " " 4n 2 (iV + 8)3 e « [( - 1) T(n + l)r(3n) + ~ ^ + 8 ) 
+4n(14iV + 40)D(n)], (114) 

a 1 3 f i ^ + 2 ) , 2 fr i^i 12 ^ r (2n) 2 o (N | o, 

P ""2 2n(iV + 8) n + 4n 2 (iV + 8) 3 nK ' r(ra + l)r(3ra) [ + J 
+n(14AT + 40)D(n)], (115) 

* -3 i e " i e ' rr.v i .^ 2 i ( ir 4n(iV + 2)r(2n)2 i riim 

5n " 3 + n + 2n 2 (iV + 8) 2[(iV + §) + I> + l)r(3n) ] ' (U6) 

Two important features emerge from this exact picture of the critical exponents shown 
above. It is instructive to discuss the exponent r) n . First, the sign of the lowest (0(e 2 )) 
contribution to the exponent r/ n is determined by the value of n. Remember that this fact 
merely reflects the change of sign of the two-point diagrams depending on the value of n. 
Second, instead of a global factor proportional to n the exact solution has a dependence on 
n coming from a product of V functions having n on their arguments as well as a finite sum 
with terms which depend on n. This property is valid for all critical exponents and appears 
explicitly at two- and higher-loop corrections. This is quite a remarkable new feature of 
arbitrary isotropic competing systems, going beyond the simple polynomial dependence on 
n found using the orthogonal approximation. 

The universality class reduction to that from the Ising-like one in the limit n — > 1 is 
obvious. Moreover, the case n = 2 correctly reproduces the solution of an earlier two-loop 
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calculation. Actually, the results given above extend the calculation of rj n for arbitrary 
n to include 0(e 3 ) corrections. In addition, all critical exponents presented here at least 
up to two-loop order generalize all previous results for arbitrary higher character isotropic 
Lifshitz points. Therefore, at the loop order considered the results above represent the 
complete solution to the critical exponents of the CECI model for arbitrary types of isotropic 
competing interactions. 

A numerical analysis for comparison between the results obtained either using the or- 
thogonal approximation or in the exact calculation are worthwhile. The case n — 2 will be 
analysed later. Here we display the numerical results for the cases n — 3, 4, 5, 6. 

It would be appropriate to calculate the exponents in either approach in a particular 
case in order to see if the difference is meaningful. The usual e-expansion is good enough 
for the numerical estimation of critical exponents associated to three-dimensional critical 
systems even though the expansion parameter is not small. We can then ask ouselves if the 
same analogy is valid in order to extract concrete results from a fixed value of the space 
dimension and number of components of the order parameter field for arbitrary isotropic 
higher character Lifshitz points. We shall look at space dimensions which yields e n = 1 in 
analogy to the method used to extract numerical results from the ordinary e— expansion for 
three-dimensional systems. 

For an eleven-dimensional lattice, take N — 1 and n = 3 for the isotropic third character 
behavior. The orthogonal approximation yields rj 3 = 0.006, v 3 = 0.177, a 3 = 0.046, (3 3 = 
0.446, 73 = 1.064 and 5 3 = 3.385. The exact calculation produces the exponents 773 = 0.002, 
v 3 = 0.174, «3 = 0.085, (3 3 = 0.435, 73 = 1.046 and $3 = 3.390. The maximal deviation 
(4.1%) occurs for the exponent a 3 followed by the deviation in the 73 exponent (1.8%) and 
an error in the exponent (3 3 around 1%. The other exponents have deviations smaller than 
0.5%. 

Consider the case N = 1, d = 15, n = 4. The results from the orthogonal approximation 
are: 774 = 0.005, v 4 = 0.131, a 4 = 0.036, /3 4 = 0.459, 74 = 1.046 and <5 4 = 3.279. The 
exact calculation yields 774 = -0.001, v 4 = 0.129, a 4 = 0.058, /? 4 = 0.449, 74 = 1.029 and 
64 = 3.282. The maximal deviations takes place in a 4 (2.2%), 74(1. 7%) and /? 4 (1%), while 
the other exponents have errors smaller than 0.6%. 

Let us examine the case N = 1, d = 19 and n — 5. The isotropic exponents in the 
orthogonal approximation are 775 = 0.004, v§ = 0.104, a§ = 0.030, ^5 = 0.467, 75 = 1.036 
and S 5 = 3.219. The exact exponents are 775 = 0.001, z/ 5 = 0.102, a 5 = 0.064, (3 5 = 0.475, 
75 = 1.020 and 5 5 = 3.220. The maximal deviations are in a 5 (3.4%), and 75(1.6%), whereas 
the remaining exponents have errors smaller than 0.8% 3 . 

This analysis leads us to conclude that the orthogonal approximation is very precise 
to predict numerical values in specific situations, since the deviations are negligible when 
compared with the exact calculation. Moreover, the above data indicate that the deviations 
are under control no matter how the number of neighbors increases. 

The extra insight from the field-theoretic viewpoint is that the more neighbors are intro- 
duced and coupled trough isotropic competing interactions the more the space dimensions 



3 The maximal deviations also occur for a n and j n for n = 6, 7, 8 and decreases for increasing n. 
For n = 8, the maximal error for o§ is about 2.7%, whereas the maximal deviation for 73 is 1.3%. 
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seem to split. Then, one line with competing interactions up to second neighbors behaves 
for all practical purposes as having 2 dimensions. Pushing the argument further, a line with 
n neighbors interacting through competing interactions seem to have n dimensions. This 
is a striking general property of the field theory under consideration: in the massless limit 
presented here, when the free critical propagator is proportional to a 2nth power of momenta 
each space direction "splits" in n dimensions. This is a kind of degeneracy in the dimen- 
sion, which can only be unveiled when more participants (neighbors) are allowed to take 
place in the isotropic competing system. This is a new aspect of systems whose Lagrangians 
have kinetic terms described by higher derivatives of the field. Further implications of this 
phenomenon remain to be investigated. 

Now, let us show that these findings easily reproduce and extend the original calculation 
done by Hornreich, Luban and Shtrikman [3] for the n = 2 case. 



VIII. EXACT ISOTROPIC EXPONENTS FOR THE SECOND CHARACTER 

CASE 

In the last section we found the critical exponents for the isotropic case when competing 
interactions among arbitrary neighbors are allowed. We now discuss its reduction for the 
usual second character Lifshitz critical behavior. Since we want to compare our results 
with other which already appeared in the literature we shall derive the critical indices using 
normalization conditions and minimal subtraction of poles. 



A. Critical exponents in normalization conditions 

The fixed point can be calculated by replacing n = 2 in Eq.(107) in order to obtain: 

6e 2 . 1 r (41iV + 202) 1.. . , 

The renormalization functions 74,(2) (1*2) and 74,2(2) (1*2) can be expressed in the following 
simple manner in terms of u 2 : 

(iV + 2) r 131 . 9 29(iV + 2)(iV + 8) „ , , 

lm)M = -- — — 1 + e 2 ii! + — -uh 118 

w n ; 240 L 120 J 2 28800 2 v ; 

%H2)(u 2 ) = u 2 [l - ^e 2 + -^u 2 \. (119) 

Substitution of the fixed point in the first equation gives directly the anomalous dimensions 
772. Thus, we get 

3(N + 2) (iV + 2) (41AT + 202) 23 

12 20(AT+8) 2 2 10(iV + 8) 2[ 10(iV + 8) 2 80 J 2 ' 1 ; 

From the scaling law relating r] 2 with v 2 we obtain: 



31 



1 (N + 2) (iV + 2)(15iV 2 + 89iV + 4) 

^ 2 " 4 + 16(iV + 8) e2 + 960(iV + 8)3 62 (121) 

Using the scaling relations to obtain the remaining critical exponents, we find: 

(AT + 2) (iV + 2)(15iV 2 + 98iV + 76) 2 , 

72 = 1 + 4(JV>8) 62 + 240(iV + 8)3 £2 ' (122) 

(4-AQ , (iV + 2)(-15iV 2 + 62Ar + 952) 2 

" 2 " 4(iV + 8) e2 + 240(iV + 8)3 c * ; (123) 

_ 1 3 (iV + 2)(80iV + 514) 

^ 2 " 2 " 4(iVT8) 2 240(iV + 8)3 621 {1M} 

1 (5iV 2 + 867V + 332) 

52 = 3 + 2 62 + 40(iV + 8) 2 £2 - (125) 

Since the earlier calculations in the literature were performed using minimal subtraction, we 
shall discuss it next. 



B. Critical exponents in minimal subtraction 

Now, let us obtain some renormalization functions at the fixed point using minimal 
subtraction. Since these objects are universal, finding them is equivalent to find the fixed 
point. Requiring minimal subtraction of the renormalized vertex r^L, the cancellations 
among logarithmic integrals for arbitrary external momenta indeed take place as expected. 
It leads to the expression of the bare dimensionless coupling constant u Q2 written in terms 
of the renormalized dimensionless coupling constant u 2 , namely 

r „ (TV + 8) AN + 8) 2 (41N + 202), ,. 

From the definitions of the normalization constants, we can write the following expressions 

(iV + 2) 2 

%0 = 1 + ~^^ U 2 



2 



(jV + 2)(7V + 8) 1 23 
+ 4320 [ e 2 518400e 2 J 2 ' 10 
(N + 2) 

V(2) = 1 + U 2 



(N + 2) (AT + 5) AT + 2 
+ i 36e 2 + 144e 2 J 2 ' 



Consequently, the functions 74,(2) (^2) and 7^2(2) (^2) read 
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. (N + 2) 2 23(iV + 2)(iV + 8) 3 

7<a^(«2) = -« H (129) 

W2n ; 240 2 172800 2 ' V ; 



V(2)(«2) = ^^«2[1 + ^tia]. (130) 

The fixed point can be calculated and shown to be 

, 6 (4UV + 202) 2 

M2 = (jvTsf 2 " 5(iV + 8)3 ^ (131) 

Replacement of this expression for the function 74,(2) (^2) gives precisely the exponent rj 2 of 
the last section up to O(eij). The corresponding expression for 74,2(2) (t^) is given by 

-* (iV + 2) f (iV + 2)(13iV + 41) 

The last expression is actually the same as that coming from normalization conditions, 
therefore leading to the same exponent v 2 whereas the remaining exponents are obtained 
using the scaling laws. Thus, the equivalence between the two renormalization schemes is 
complete. 



C. Discussion 



First of all, our results for the isotropic n = 2 case generalizes those in the seminal 
work by Hornreich, Luban and Shtrikman [3]. To see this, we make the identifications 
e a = €2 , ?7z4 = r]2 and 1^4 = z/ 2 . The equations (10a,b) in [3] are identical to our results for 
7)2 and V2 Eqs.(120)-(121). The step forward within our method is the exponent rj 2 which 
is obtained up to 0(e 2 ) for the first time. Furthermore, using the scaling relations reported 
in our previous letter [15], we found all critical exponents (Eqs.(122)-(125)) exactly, at least 
up to two-loop order which constitutes another new result for the usual isotropic case. 

Next, let us confront the results coming from the generalized orthogonal approximation 
with those from the exact solution at the same loop order. Both the exact and the approx- 
imated one-loop exponents are the same. This can be seen from the n = 2 particular case 
or from the generic n isotropic criticality by the direct examination of the explicit results 
shown in the present article. Two-loop deviations start at n = 2 and for higher n. 

Consider the case n — 2. Take a magnetic system which has N = 1 in a seven-dimensional 
lattice and analyse the exponents in each case. We shall restrict ourselves to three significa- 
tive algarisms in our discussion. First use the orthogonal approximation. At two-loop order, 
the correlation length exponent yields the result u 2 = 0.276, and the anomalous dimension 
is given by r] 2 = 0.009. The susceptibility, specific heat, magnetization and magnetic field 
exponents are given by 7 2 = 1.103, a 2 = 0.061, (3 2 = 0.418, 5 2 = 3.616. Now use the 
exact two-loop calculation. We obtain the following numerical values for the critical indices: 
v 2 = 0.271, 7/2 = -0.006, 72 = 1.087, a 2 = 0.100, (3 2 = 0.406, S 2 = 3.631. 

Therefore, for all the exponents the difference using either approach starts in the second 
significative algarisms. Specifically, the maximal error made by using the orthogonal approx- 
imation takes place for the specific heat exponent with a deviation of 3.9% , whereas the 
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minimal error occurs in the correlation length exponent whose difference is approximately 
0.5% . This is a strong evidence that the orthogonal approximation is very good to give 
reliable information for the isotropic case in this specific situation. This feature was already 
encountered for uniaxial anisotropic cases, where the approximation showed its usefulness 
for three-dimensional uniaxial systems. We hope that these results may stimulate the search 
for these exponents using Monte Carlo numerical simulations in the particular case of second 
character isotropic Lifshitz critical points. 

The remarkable agreement between the numerical results for the critical indices above 
mentioned either using the orthogonal approximation or the exact two-loop calculations 
corroborates previous conjectures that the orthogonal approximation is not only good to 
describe uniaxial systems pertaining to the anisotropic second character Lifshitz critical 
behaviors but also arbitrary isotropic higher character Lifshitz points. 

In fact the numerical analysis can be carried out for higher-dimensional lattices for higher 
values of n. Extending this argument, the case L = An, N = 1 and d = L— 1 yields the same 
value for the expansion parameter and should not deviate very much when both calculations 
are compared. Perhaps the study of the most general arbitrary isotropic points via numerical 
tools might be worthwhile as well, since now we have satisfactory numerical results coming 
from a purely analytical field theoretical investigation. 

IX. CONCLUSIONS AND PERSPECTIVES 

In this paper we have discussed the field theoretic description of the most general compet- 
ing system, which has a simple lattice model representation named CECI model. It consists 
of a modified Ising model presenting the most general type of competing exchange couplings 
among arbitrary neighbors and includes other models previously reported. We have derived 
explicitly the critical exponents in the anisotropic as well as in the isotropic situations at 
least up to two-loop level. The CECI model and its field theory representation generalize 
the description of the second character Lifshitz universality classes in a nontrivial way. In 
particular, strong anisotropic scale invariance is exact up to the loop order considered here. 
The universality class reduction is a general property of both anisotropic and isotropic crit- 
ical behaviors. It implies that when the interactions beyond first neighbors are turned off, 
the Ising-like universality classes are recovered. This feature is manifest in all exponents. 

The anisotropic exponents were calculated by using the generalized orthogonal approx- 
imation. The calculation of loop integrals is consistent, since it is rooted in the physical 
property of homogeneity. It is required for a satisfactory renormalization group treatment 
with several independent relevant length scales represented by each correlation length associ- 
ated to the several competing axes. The fixed point is the same irrespective of the competing 
axes under consideration. In close analogy to the second character case, this result is ex- 
pected to hold in arbitrary perturbative orders. The second character Lifshitz exponents 
are easily recovered as a particular case of the generic anisotropic situation described in 
the paper. Although it is desirable to have an approach that does not require the use of 
approximations for the anisotropic behaviors, in the present moment it is not obvious. It is 
a consequence of the appearance of many competing subspaces simultaneously which makes 
the exact calculation (if not impossible) very difficult. We hope the techniques developed in 
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the present work shed light on a quest for a solution without the necessity of approximations 
for the anisotropic cases. 

The isotropic behaviors were calculated using two different trends to evaluate the Feyn- 
man diagrams. The first of them makes use of the orthogonal approximation. It can be 
noticed that the isotropic behavior cannot be obtained from any type of anisotropic behav- 
ior within the framework of the e^-expansions developed in the present work even though 
the same kind of approximations are employed in both cases. This generalizes the previous 
situation taking place for second character Lifshitz points. Next, we attacked the diagrams 
without making any approximation. The result is that deviations in the calculation of the 
critical exponents start at two-loop level in comparison to the outcome provenient of the 
orthogonal approximation. We obtain as a particular case of the generic isotropic behavior 
the second character isotropic behavior, which extend the results first derived by Hornre- 
ich, Luban and Shtrikman [3]. In this way, we obtain rj 2 up to three- loop order and i/ 2 at 
two-loops as well as the remaining exponents via the scaling relations derived previously in 
[15], a result which was lacking since the discovery of the second character isotropic Lifshitz 
point. 

The most immediate application of the formalism just presented is the calculation of 
all universal amplitude ratios of certain quantities close to generic higher character Lifshitz 
points. It is amazing that a detailed scale theory for these quantities is lacking in the 
literature, even for the simple second character Lifshitz points. In fact, some results were 
presented for the susceptibility [30] and specific heat [31] amplitude ratios at one- loop order. 
The latter amplitude ratio proved to be reliable to explain the experimental result associated 
to the magnetic material MnP [32] . Nevertheless, a thorough renormalization group analysis 
is necessary in order to have a better comprehension of the several scale transformations in 
each competing subspace and the role played by them in the treatment of generic amplitude 
ratios. In principle we can extend the formalism presented here for the calculations of 
amplitude ratios including the most general competing system. 

The treatment of finite-size effects for the most general Lifshitz critical behavior can be 
developed as a direct extension of the approach to the Ising-like critical behavior [33,34]. 
The applications might include systems which are finite (or semi-infinite) along one (or 
several) of their dimensions, but of infinite extent in the remaining directions. It would be 
interesting to see how different competing axes alter the approach to the bulk criticalities, 
for example in parallel plate geometries. The utilization of different boundary conditions in a 
layered geometry would be particularly simple and instructive to see how the generalization 
works for arbitrary competing systems. It could be used to investigate how amplitude ratios 
change with different boundary conditions with respect to the situation occurring in bulk 
systems [35]. 

Typical examples are systems which are finite in all directions, such as a (hyper) cube 
of size L, and systems which are of infinite size in d! — d — 1 dimensions but are either of 
finite thickness L along the remaining direction (d-dimensional layered geometry) or of a 
semi-infinite extension. The presence of geometrical restrictions on the domain of systems 
also requires the introduction of boundary conditions (periodic, antiperiodic, Dirichlet and 
Neumann) satisfied by the order parameter on the surfaces. In particular, the validity 
limits of the e^-expansion for these systems and the approach to bulk criticality in a layered 
geometry can be studied [35]. 
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One interesting aspect of the generalized orthogonal approximation is that it can actually 
address the problem of calculating Feynman integrals originating in field theories in the 
massless limit with odd (greater than 2) powers of momentum in the propagator as well. 
This subject goes beyond the realm of critical phenomena in competing systems. It might 
be useful for treating perturbatively a recent proposal of an effective quantum field theory 
with cubic kinetic terms [36]. In addition, the general framework can be used to treat 
perturbatively other effective quantum field theories with higher derivative kinetic terms 
which break Lorentz invariance in the infrared regime as the effect of combined gravitational 
attraction and repulsion [19]. This type of theory resembles very much the second character 
Lifshitz critical behavior above discussed. One can expect that introducing more and more 
competing gravitational interactions in the infrared (long distance) limit higher powers will 
appear in the kinetic terms of the corresponding effective field theory. The perturbative 
analysis of this system would be quite analogous to the generic higher character discussed 
in the present work. 

To summarize, we have described the generic higher character Lifshitz critical behaviors. 
New field theoretical tools were exposed resulting in new analytical expressions for all the 
critical indices in the isotropic as well as in the anisotropic cases at least at 0(e|). Other 
aspects like crossover phenomena and tricritical behavior for this model remain to be studied. 
New anisotropic behaviors in the absence of a uniform ordered phase for the CECI model 
are under current investigation. 
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APPENDIX A: FEYNMAN INTEGRALS FOR ANISOTROPIC BEHAVIORS 

In general the critical exponents and other universal ammounts are independent of the 
renormalization group scheme employed. The explicit calculations in this section are pre- 
sented in such a way that the results can be checked using more than one renormalization 
procedure. We shall now describe the generalized orthogonal approximation (GOA) for the 
integrals appearing in the anisotropic cases. 

In order to accomplish the task of calculating the critical indices at least up to two- 
loop level, we need a minimal set of Feynman diagrams to work with. The one- two- and 
three-loop integrals we shall need to determine are given by the following expressions 

d d ~^ m -qU L n=2 d m -k {n) 

KE!U(*<n) + K[ n) rr + (? + p?\ (EL.(*f„)) B + ? 2 ) 
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^-EL 2 ^ gi ^EL 2 ^ g2 n^ =2( i m «A; 1(n) n^ =2 ^A; 2(n) 



(ft 2 + Ei# 2 w )") (^i + Et^ 



h = 



X [(ft + ft + P) 2 + (EL,(*Kn) + *2 ( „) + K[ n) y) n } ' (A2) 

(ft 2 + Et^w^) (cp - ft) 2 + Et 2 ((^;„) - ^(n)) 2 ) n ) 

x- , , (A3) 

[qI + Et#2(n))") [(ft - ft +P3) 2 + E!U((*K») - *2(n) + *£ (n) )T] 



d d ~^n=2 m nqid d-zL2 ^ q2d d-EL 2 mn q 3 U^ =2 d mn k 1(n) 

ft 2 + Et# 2 w )") (ft 2 + Et# 2 % } ) n ) (ft 2 + Et 2 (%))") 

n^ =2 rf m "A; 2(n) n^ =2 ^A;3 (n ) 



=2 

1 



[(ft + ft - P) 2 + E„= 2 (( fc l(n) + k 2(n) ~ k[ n) ) 2 ) 



[(ft + ft - P) 2 + E„= 2 (( fc l(n) + fc 3(n) - k[ n) ) 2 ) n 



(A4) 



We stress that in the above expressions P, p 3 and p are external momenta perpendicular to 
the various competing axes, whereas K',, , kL, and k',s are external momenta characterizing 
the nth competing axes, respectively. For arbitrary values of the external momenta, these 
integrals can be calculated by making use of approximations very similar to those first 
developed to describe second character Lifshitz points [8]. As a matter of fact, the generalized 
dissipative approximation was formerly used to compute the critical exponents out of the 
anomalous dimension and correlation length exponents corresponding to space directions 
perpendicular to the competing axes for this model. Indeed, the generalized orthogonal 
approximation was figured out using a similar analogy. Nevertheless, since the dissipative 
approximation cannot approach the isotropic cases, we shall not describe it in this paper. 
Instead, we shall make use of the generalized orthogonal approximation, for it can approach 
both anisotropic and isotropic behaviors. 

Certain useful identities will be derived in order to solve the integrals above. Let us 
proceed to find out them. Our starting point is the integral derived in [8], namely 

/°° 1 771 -m 

dx 1 ...dx m exp(-a(x\ + ... + x 2 m ) n ) = —T(—)a^S m . (A5) 
-oo ^ n ^ n 

After choosing r 2 = x\ + ... + x 2 m , one can take y = r n to write last equation in the form 
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dyy™ exp{-ay ) = -o'» r (^> ( A6 ) 



The integral 
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/oo 
exp(-ax 2k - b(x + x ) 2k )dx (A7) 
-oo 

cannot be solved exactly for all k. In fact for generic k > 2 it has no elementary primitive. 
Nevertheless, one can select only the homogeneous function of o by using the following ap- 
proximation. First, make the change of variables y = x k . Second perform the approximation 
(x + x ) k = x k + Xq. Thus, this integral becomes: 

exp(-ax 2k - b(x + x ) 2k )dx = exp(-byl)- \ exp(-(a + b)y 2 - 2byy )y^ 1 dy (A8) 
oo * Jo 

Next, perform another change of variables, namely, y' = y + We then obtain: 

JZc exp(-ax 2k - b(x + x ) 2k )dx = exp{- ^y 2 )^ exp(-(a + b)y' 2 ){y' - 

- exp(-(a + b)y> 2 )(y> - ^)^dy']. (A9) 



Since the integrals are convergent, we can make use of the approximation (y' — j-)* 1 = 
y'^ 1 + .... The elipsis stands for the remaining terms that will be subtracted from the last 
integral, producing a type of error function. The original integral is then approximated by 
its leading contribution, i.e. 



J — c 



exp(-ax 2k - b(x + x ) 2k )dx = exp( °\x 2k )\v(^-)(a + (A10) 



It can be easily shown that the generalization for the m-sphere yields 

JZc exp[-a(x 2 + ... + x 2 J k - b(( Xl + x 0l ) 2 + ... + (x m + x Qm ) 2 )] k dx l ...dx m 2* (All) 

ab ~,2k\ 1 c p/ m ' 
' a+b X J2k° ml \2k' 



^(-^)sS m r(£)(a + &)"*. 



We found appropriate to name this approximation the generalized orthogonal approximation, 
for it is a natural generalization of that discussed in the usual second character case [8]. 

Let us now start our calculation of loop integrals given by Eqs.(Al)-(A4). We can 
calculate the one-loop integral using two Schwinger parameters. Using the formula derived 
above, the integration over quadratic momenta can be shown to be given by 



/2 = 2V^™/( f*— U I da ' da2 ^-^T^ 

x(a 1 + a 2 )-( d - zE ^ L ) j (nt 2 rf m "A; (n) )exp(- ai ^(^ n) r-a 2 ^((% + K; n) ) 2 r). (A12) 

J n=2 n=2 



We can now expand the argument of the last exponentials. Notice that now we have a prod- 
uct of independent integrals corresponding to the momentum components along arbitrary 
competing subspaces. Those integrals cannot be performed analytically. If we use the ho- 
mogeneity property of the remaining integrals in the arbitrary competing external momenta 
scales, we can simplify the calculation by utilizing the generalized orthogonal approximation. 
In fact, taking (k^ + K'^) n = k r ^ + K'fy and using (All), we obtain: 
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jd m -k {n) exp(- ai k^ - a 2 (k {n) + K'f n ) = ±S mn 



"(„) "2 \^(n) -r v (n)y 
1 _ (") \J 



x ( ai + a 2 )-^exp(- a -^)T(^) (A13) 



Therefore, we can write the integral as 



r _ !<; rf (d-&^) V nL ^r(^) 

^-jteJi 2 )(II "= 2 Vo X 



x exp( ^ ) («! + a 2 ) 2 . (A14) 



Take x = «i(P 2 + Et 2 ((^H) 2 )") and f = ^+EL(W 2 ) B )- Then > define « = if*- 
Consequently, the parametric integrals can be performed with this change of variables. Using 

the identity 

T(a + bx) =T(a) [ 1 + bx ip(a) + 0(x 2 ) 1 , (A15) 
and expressing everything in terms of the parameter leads to the following expression for 

TV m n ' 



7 * = i[VE^™ n) r(2 - Et 2 ^)(nL 2 ^£^)](i - W - Et 2 trMf ) 

x £ d«( v (i - v)(p* + Et 2 ((^ n) ) 2 ) n ))^- (A16) 

This is a homogeneous function (with the same homogeneity degree) in (P,K^) . The 

factor [S (d _ E L_ 2mn) T(2 - E„ =2 ff)( n «=2 2n )1 is S oin S to be absorbed in a redefinition 
of the coupling constant after performing each loop integral and shall be omited henceforth. 
We can follow two routes from last equation. The first one is to perform the v integral in 
terms of products of the Euler T functions. It will be useful in the calculation of higher- 
loop integrals since we need the subdiagrams of this type in order to compute the complete 
integral. Then, we obtain 

HP,K[ n) ) = (P* + Yi«K{ n) m)^± (l + h mL e L ) , (A17) 

n=2 L 

This is appropriate to calculate the critical exponents only using normalization conditions. 
But we would like the solution in a form suitable for minimal subtraction as well. The 
second possibility convenient for both types of renormalization schemes is to expand the 
last integral as 

f dv( V (i - v)(p 2 + ^((K'^rr))^ = 1 - € -§m n»)> ( A18 ) 

J ^ o 



n=2 



where 



L(P,K{ n) ) = fdv ln[v(l-v)(P" + J2((K{ n) ) 2 y% (A19) 

J n =2 
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Hence, this integral reads 

I 2 (P, K{ n) ) = 1 (l + (h mL - l)e L - e j-L(P, K{ n) )) , (A20) 

where h mL = 1 H 2 "~ 2 2ra — . Notice that whenever m 3 = ... = mi = 0, /i m2 = 

[i2]mj and the usual anisotropic Lifshitz critical behavior is trivially obtained from this more 
general competing situation. This form is convenient for the renormalization using minimal 
subtraction. Instead, for normalization conditions we have: 

hsp, = ■■■ = hsp L = — ( 1 + h mL e L ) , (A21) 

since L{SP 1 = ... = SP L ) = -2, with SP 1 = (P 2 = 1, K[ n) = 0),..., SP L = (P = 0, (K ( ' L) ) 2 = 
!)• 

The simplifying condition (fc( n ) + K'i n \) n — + -^(™) for the one-loop integral can be 
generalized to the higher-loop graphs. It is translated in the statement that the loop momenta 
characterizing a certain competition subspace in a given bubble (subdiagram) do not mix to 
all loop momenta not belonging to that bubble. The simplest practical application of this 
principle can be viewed in the calculation of the "sunset" two-loop integral I 3 contributing 
to the two-point function 



(tf + EL»(*?(„)) B ) (? 2 2 + £t 2 (%) 



x f ; , (A22) 

[(9i + Q2 + P? + (S = #i W + h {n) + K' [n) m 

Defining K'!-, = k^ n ) + K',-, and using the condition k 2 ^ n yK'/-, = 0, one can solve the 
integral over q 2 , k 2 i n ) first, picking out only the homogeneous part of each individual integral. 
The remaining parametric integrals contains the divergence (pole in e£) and can be solved 
as before. Using Eq.(A17), we obtain: 

I 3 (P,K') = -(l + h mL ) 

x I d d -^ m - qi ui =2 d m -k l{n) 

J (q! + Et#W") [(ft + p) 2 + £L,((*Kn) + K{ n) ) 2 ) n P ' 

Using Feynman parameters, integrating the loop momenta along with the remaining para- 
metric integrals, and expanding the resulting Y functions in we find: 

I 3 (P, K') = (P 2 + jr < 2 ")^(l + 2h mL e L - \e L - 2e L L 3 (P, K{ n) )\ (A24) 

n=2 

where 



L 3 (P, K') = f dx(l - x)ln[(P 2 + V K' n 2n )x(l - x)]. 

J ° n=2 



(A25) 
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At the symmetry points SP n , it can be rewritten as 



I 3SPl = ... = h SPL = — (1 + 2h mL e L + |e L ). (A26) 



W= ^) = - = W= ^r) = 8^(1 + 2hm L eL + ~ A e L ). (A27) 



From the above equation we can derive the expressions: 

To complete our description of the 1PI two-point vertex parts, consider the integral 

d d -^n=2 m nqid d-^ =2 m nq2d d-Z^ =2 mn q 3 U^ =2 d mn k 1{n) 

(ft 2 + Et 2 (%))") (<Z2 + E;= 2 (*l(n)) n ) (ft 2 + Et#3 2 (n)) n ) 

n^ =2 (i m "A; 2(n) n^ =2( i m «A;3(„) 

X 



[(ft + ft ~P) 2 + ££= 2 (( fc l(n) + k 2 (n) ~ k') 2 )' 



x = — , (A28) 

[(ft + ft -P? + EL((*1W + fc 3(n) - fc( B) ) 2 ) B ] 

which is the three-loop diagram contributing to the two-point vertex function. Incidentally, 

there is a symmetry in the dummy loop momenta q 2 —> ft and k 2 ( n ) — > &3( n ). Concerning the 

integrations either over q 2 , k 2 ( n ) or q 3 , k^^, we use the condition (&2( n ) + (&i — -^')) n = ^2(n) + 

(^i(n) — -^( n )) n wnen the integration is performed over k 2 as well as (fc 3 ( n ) + (fci( n ) — -^( n ))) n = 

^3(n) + (^i ~~ ^'(n)) n wnen the integral over k 3 is realized. The two internal bubbles, which 

are represented by the integrals over (ft,/c 2 ( n )) and (ft,&3( n )), respectively, are actually the 

same, resulting in i 2 ((ft — P), (&i(n) — ^L)))- Next ta ^ e ^* ~~ * ~~ ^> ^(n) ~~ * ~^[n)- Using a 
Feynman parameter and proceeding in close analogy to the calculation of -Z3 we find: 

L -1 

J B (P, ffJ B) ) = (F 2 + ^ < 2n )^(l + 3^e L - c L - 3e L L s (P, K{ n) )l (A29) 



n= 2 - L 



At the symmetry points SP\, SPl one obtains: 



^) = ... = W= §S) = J(i + 3^ ei + i C£ ). 

(L) L 



Finally we compute one of the two-loop diagrams contributing to the four-point function, 
namely 

d d -^n=2 m nqid d-E^ =2 ™-g 2 n^ =2 (i m " k 1(n) U^ =2 d m " k 2{n) 



(ft 2 + £t# 2 w ) n ) ((^ - ?i) 2 + Et 2 ((^„) - *1<»))T) 

x- , . (A31) 

ft 2 + Et# 2 % } n [(ft - ft + p 3 ) 2 + Et 2 ((^(n) - h(n) + K (n) m 



Recall that P = pi+p 2 , p% {% = 1, 3) are external momenta perpendicular to the competing 
axes. On the other hand, K'^ = k[^+k 2 ^, and k' i( ^ (i — 1, 3) are the external momenta 
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along arbitrary competition directions. We can integrate first over the bubble (#2, &2(n))- It 
is convenient to choose Schwinger parameters in the calculation. Then, we use two Feynman 
parameters and solve for the loop momenta to obtain the following parametric form 

■I "> 7 — \ (•!„ 2 - 



2'"" ^ r (w \ ; ' J 2/; 



: f dyy(l-y)^~ l f dz\yz(l - yz)(P 2 + £ ffg)) + y(l - y){p\ + £ fc^) 

^0 ^ n=2 n=2 

L 

2yz(l - y)(ps.P + £(-!)"*;?„) J#))l " (A32) 



n=2 



The integral over y is singular at y — 1 when = 0. Replace the value y — 1 inside the 
integral over z [1,8], integrate over y and expand the Gamma functions in e^. This implies 
that 

h = 2? ( X + 2 " ^ L ~~ 6lL(jP ' ' (A33) 

This form is particularly suitable for the renormalization procedure using minimal subtrac- 
tion. For the purpose of normalization conditions, the value of this integral at the symmetry 
points discussed before is given by 

I ASPl = ... = h SPh = (l + 2 h mL e L + ie L ) . (A34) 

The method proposed here is equivalent to a new regularization procedure to calculate Feyn- 
man integrals whose propagators have any combination of even powers of momenta. We can 
define the measure of the m„-dimensional sphere in terms of a half integer measure. In 
fact, taking k = p 2n , one has d rrln k = d^p = -^jf^~ x dpd^l mn . Hence, the approximation 



required to solve the integrals results that the new "measure" d^p is invariant under trans- 
lations p' = p + a. This is a simple generalization of the same property valid for the usual 
m 2 -fold Lifshitz behaviors. 



APPENDIX B: ISOTROPIC DIAGRAMS IN THE GENERALIZED 
ORTHOGONAL APPROXIMATION 

The computation of the Feynman integrals using the generalized orthogonal approxima- 
tion is simpler in the isotropic cases, since there is only one subspace to be integrated over. 
At the Lifshitz point 5o n = r nn < = and solely the 2Lth power of momentum appears in the 
propagator for the case of Lth character isotropic critical point. The isotropic analogous of 
the one-loop integral contributing to the four-point vertex part is 

/ dm " k 
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We can use two Schwinger parameters and the orthogonality condition (k + K') n = k n + K' n , 
resulting in the expression 

/poo poo 
d m "k / / da 1 da 2 e-( ai ^X*re- 2a * K ' nkn e- a « K ' a r. (B2) 
Jo Jo 

Turning to polar spherical coordinates, take r 2 = k\ + ... + k 2 mn . Making the transformation 
k n = p the volume element becomes d mn k = ^p^^dpdQ^ = d^p. The former integral 
with a nth power of momenta changes to a quadratic integral over p. After discarding the 
infinite terms which change the measure d~^k under the translation y' = y + only the 
leading contribution is picked out and we have 

/ d m -ke- a ^ n - bkn = [ d^pe-*-* = a-^e^ir(^)5 m „. (B3) 

J J Zn 

When this result is replaced into the expression of I 2 (K'), we get to 

J 2 (K') = ^[1-^(1 + ^' 2 ))]. (B4) 

Henceforth we absorb the factor of S rrin in this integral through a redefinition of the coupling 
constant and shall do so after performing each loop integral for arbitrary vertex parts. Note 
that this absorption factor is different from that arising in the anisotropic case in the limit 
d — > m n = An — e n . Since the geometric angular factor coming from the anisotropic cases 
becomes singular in the above isotropic limit the attempt of extrapolating from one case 
to another is not valid, at least within the framework of the ex-expansion presented in this 
work. This is a further technical evidence that the isotropic and anisotropic cases have to 
be tackled differently. Thus, 

J 2 (/0 = ^[1-^(1 + (B5) 

The suitable symmetry point (K' 2 ) n = 1 useful for the purpose of normalization conditions 
leads to the following simple outcome 

(B6) 



The next step is the evaluation of the integral 

d mn hd mn k 2 



((h + k 2 + K'f) n {kiY{kiy 



(B7) 



Integrate first over k 2 . Take K" = k l + K' and use the condition (k 2 + K") n = k% + K" n to 
obtain: 

h = h l+ ^ ] J m+Unmr- (B8) 

Utilizing a Feynman parameter, we can integrate over k\. After the expansion m n = An — e n 
is done inside the argument of the resulting T functions and using the expression 
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d^q „ 1 r(^)r(a - ^)(m 2 - A; 2 )^-^ m „ 



(q 2 + 2k.q + m 2 ) a 2n T(a) 
the integral J 3 can be found to be 

(K' 2 ) n 1 2 



(B9) 



8ne„ An n 

At the symmetry point, this reduces to 



[l + e n (- ^ 3 (^' 2 ))]. (BIO) 



implying that 



^ l + h^ (B11) 



! ' flP = / s = -- + (B12) 



d(K l2 ) nX ' 3 6 8ne 4 „ L 4n 
The 3-loop integral J 5 is given by 



((fci + A; 2 + tf')T((*i + ^3 + ') 2 ) B (*?) B (*2) B (*§) B 



(B13) 



where we took for convenience the redefinition K' — > — X'. The integrals over k 2 and A; 3 
are the same. Thus, following analogous steps and employing the same reasoning as in the 
calculation of J 3 we get to 

(K' 2 V l 1 3 

' l " [l + e n (^--L 3 (K"))}. (B14) 



6ne 2 2n n 

At the symmetry point, the following expression follows trivially 

dh , 1 r- 2 



The two-loop integral Z4 in the isotropic behavior is 

where K' = k[ + A;f>. The integration can be done along the same lines of the computation 
performed for its anisotropic counterpart. It is straightforward to show that 

h(K") = ^[l-^-(l + 2L(K' 2 ))]. (B17) 

At the symmetry point the integral can be rewritten in the form 

, 4 (A- = l) = ^[l + § ? . (B18, 

As was shown above, these results are a natural generalization of those originally developed 
for the second character Lifshitz points. It can be checked that all integrals reduce to the 
usual (f) A values for n — 1 and reproduce the results from [8] in case n — 2. 



= / 5 = -7- T [l + -*>]• (B15) 
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APPENDIX C: ISOTROPIC INTEGRALS IN THE EXACT CALCULATION 



An interesting feature of the isotropic case is that it can be calculated exactly. We now 
proceed to yield the exact solution to the Feynman diagrams without performing approxi- 
mations. 



/ 



d mn k 



((k + K') 2 ) n (k 2 ) n ' 

Using a Feynman parameter and making the continuation d = m n = 4n — e n we get 

T(2n-f)T(f)S mn l T(n)T(n) e 
h{K } = 2TWF) " 2 Ln(K )] ' 

where L n {K') is given by 

L n {K') = [ dxx n -\l - x) n -Hn[x{l - x)K' 2 ]. 
Jo 



(CI) 



(C2) 



(C3) 



This integral is the analogous of the integral L(K') appearing in the orthogonal approxima- 
tion. Here it depends explicitly on n, and that is the reason we have included a subscript 
in it to emphasize this dependence. The integration over x together with the e n expansion 
of the Gamma functions results in 

W = - <± mn) - m + J^UK'))]. (C4) 

As before, we absorb the factor S mn in a redefinition of the coupling constant. We need to 
do this for each loop integral. Thus, 

hm = i[i - | W(2n) - m + j^MiO)]- (C6) 

This is a useful result for doing minimal subtraction. Defining the quantity D(n) = 
\ijj{2n) — ip(n) + at the symmetry point K' 2 = 1 the integral turns out to be 

hsp = -[l + D(n)e n ]. (C6) 

Now, let us calculate the integral 1%. As before, take K" — k± + K' and solve for the 
internal bubble k 2 using Feynman parameters solving the momentum independent integrals 
over the Feynman parameters and expanding the Gamma functions in e„, we end up with 

- ^-"" vSn+v h 1 + E " (B " - (C7) 



where 



_ r(2n)r(n) 

An ~ r(3n) ' (C8) 
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2n-l 



Br, 



D(n 



E 



j 2 ^ p ^p 

p=l i p=l 

2n-l (2n-l)!(-l)P +1 
p=0 2p!(2n-l-p)!(n+p) 2 



(C9) 



L 3n {K')= [ dxx^-^l-x'r-Hnixil- x)K' 2 }. (CIO) 

Again, this integral depends explicitly on n and should be compared with its counterpart 
arising in the orthogonal approximation. We then learn that for massless propagators with 
arbitrary power of momenta, the external momentum dependent part of the Feynman inte- 
grals generalizes the standard </> 4 in the manner prescribed above. 

At the symmetry point, the integral ^3 simplifies to the following expression: 

41 [n + 1)1 {on)e n 4 n 

„ 3rt— 1 2n— 1 n— 1 

p + %£ -E )+|E^i- (cm 

p=3 1 p=2 1 p=l 1 

Notice that while the first term inside the parenthesis contributes for arbitrary values of n, 
the last factor of 0(e n ) into the brackets are corrections which contribute solely for n > 2. 
(A similar feature will also take place in the calculation of I A and I 5 .) Therefore, taking the 
derivative with respect to K' 2n at the symmetry point produces the result 

„ 3rt— 1 2n— 1 n— 1 

p=3 y p=2 F p=l y 

We now calculate the three-loop integral 1$. Proceeding analogously, we can show that 
it has the solution 

t u t2n( -|\n r 2 (2n) 1 3L 3ra (i : C / ) /pio\ 

Js = K 3 r(3n)r(n+i)4 [1 + e " (a " ^4^ )] ' (C13) 

where 

2n-l 



C n = 2D(n)-iV - + - 

v ; 2 ^ p 2^ p 

p=i p=i 



E2n-1 (2to-1)!(-1)p +1 

p=0 p!(2n-l-p)!(n+p) 2 ^ C14 ^ 

At the symmetry point this result gets simplified. Taking the derivative with respect to 
the external momenta we get to: 



46 



J '»'- = (- 1 > V( B + ( S(8»W (1 + J( " ) + 1)<l! 



3n— 1 2re— 1 n— 1 



Notice that the 0(e n ) terms inside the bracket gives a nonvanishing contribution only for 
n > 2. 

To conclude, let us calculate the integral J 4 . The integral over the bubble k 2 can be 
solved directly by taking the effective external momenta as K" = —k\ — k' 3 . It has fyK") as 
a subdiagram. Using the information obtained in calculating I 2 and working out the details 
we find the intermediate result: 

r(2n-f)r(e n ) f 1 
h = fM I (hl " A " ] 

[ dz[z{\ - z)] n - l [yz(l - yz)K' 2 + y(l - y)k' 3 2 - 2yz{\ - y)K' '.k'^ . (CI6) 
Jo 

The situation resembles that in the calculation of I4 using the orthogonal approximation. 
Again, set y — 1 in the integral over z and carry out the integral over y independently. 
Performing the integral over y, we find: 

T(2n - f )r( £ „)r(2n) ,„ ,.„_,,„ „, D ,. 2l _,„ 



'< = f^^mnMfw I dz[z{l - z)]n W1 - (C17) 

Then for the purposes of minimal subtraction, it can be expressed as 

/4 ^ [1+(D( „ ) _ 1 _E^ K _ e ;f_^ ;] . (C18 ) 

We emphasize that the last 0(e n ) term in this expression only contributes for n > 2. At the 
symmetry point, we find 

2n-2 n-1 

= ^[1 + (D(n) + l)e n + — - 2 £ — )]• (C19) 

n p=i p=i 

Once again, the last two terms of 0(e n ) containing the sums in the above expression cor- 
respond to corrections in case n > 2. Since the corrections are absent when starting from 
the scratch for the n — 1 case, by neglecting them in the above expressions it can be eas- 
ily checked that all of these integrals reduce to the values of the ordinary A0 4 in the limit 
n — > 1. We learn that a general feature of the exact calculation is that higher loop integrals 
generally receive further contributions to the subleading singularities for n > 2. 

In order to make contact with the more concrete case of the usual second character 
Lifshitz point obtained by Hornreich, Luban and Shtrikman [3], we discuss the particular 
n = 2 case next. 
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1. The n=2 case 



Let us now analyse the Feynman integrals involved in the calculation of the critical indices 
of the ordinary second character Lifshitz critical behavior. This is a mere particular case of 
the most general isotropic CECI model discussed in the previous subsection. Nevertheless, 
the discussion of this particular case is useful when comparing with the original previous 
result obtained by Hornreich, Luban and Shtrikman about three decades ago. Needless to 
say, both results agree for the exponents 77^,4 and u L4: , which in the notation of Section IV 
correspond to r\ 2 and v 2 . Moreover, our treatment permits to obtain two new results for 
this behavior: the results for r\ 2 are extended including corrections up to 0{e\) whereas the 
remaining exponents are obtained through the complete set of scaling relations derived in 
[7] up to O(el). 

Since the calculation was already indicated in the last subsection, we simply quote the 
results. For calculating the exponents using minimal subtraction the most appropriate form 
of the integrals are given by 

I 2 {K>) = -[I - ^ - 3e L L 2 (K% (C20) 
h = K ' A ^-\ l ~ + 20L 32 (K'))}, (C21) 



1 r . ,7 



h = K 20? [1 " £l( 60 + 30L32 ^' ( C22 ) 

L 

/ 4 = ^[l-( T | + 6L 12 (K' 2 )) eL ]. (C23) 

L 

For the use of normalization conditions, however, it is convenient expressing these inte- 
grals at their symmetry point. Instead of calculating ^3 and I§ at the symmetry point, we 
need their derivatives with respect to K 4 at the symmetry point. Thus, we find 

hsP = -[I ~ j,]- (C24) 

/' =-2-[i + t2i eL ] (C25) 
3SP 80e L [ 120 J V ; 



5SP 



20e 



M 26 1 



1 r 1 1 

M 11 ~ i etl 



(C26) 



(C27) 
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FIGURE CAPTIONS 



FIG. L The simplest example of the CECI model with uniaxial competing interactions between 
second neighbors as well as uniaxial couplings between third neighbors. This system has three 
independent correlation lengths and presents a generic third character anisotropic Lifshitz critical 
behavior. Note that by turning off the interactions among second neighbors leads to the simpler 
third character critical behavior. 

FIG. 2. The phase diagram of a typical uniaxial second character Lifshitz critical behavior. The 
dashed lines indicate a first order transition between the uniformly ordered and modulated ordered 
phase which terminates at the Lifshitz point of second character. The parameter p is defined by 

P=f 1 - 

FIG. 3. The superposition of the two two-dimensional independent phase diagrams for the sim- 
plest CECI model. In this two-dimensional picture, the confluence of the two distinct modulated 
phases, the ferromagnetic and the paramagnetic phases occurs at the generic third character Lif- 
shitz point. 
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